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FOREWORD 


This  report  has  been  prepared  by  the  Elkton  Division  of  Thiokol  Chemical 
Corporation  for  the  Rocket  Propulsion  Laboratory  (DGSMD)  at  Edwards,  California. 

Lt.  Harold  W.  Gale  is  the  Air  Force  Project  Officer.  The  report  is  based  on  the  work 
accomplished  during  the  period  from  January  1,  1963,  to  April  1,  1963,  under 
Supplemental  Agreement  No.  6  to  Contract  AF  33(616)-6530,  an  engineering  research 
program  to  develop  advanced  methods  for  improving  the  state  of  the  art  of  solid  propellant 
grain  design. 

The  principal  contributors  to  this  report  are  Messrs.  W.  G.  Andrews,  D.  H. 
Frederick,  F.  E.  Moore,  D.  Saylak,  A.  Stornelli,  D.  D.  Thomas  and  R.  H.  Thompson. 
Mr.  James  F.  Hoebel  is  the  Program  Manager. 
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ABSTRACT 


The  work  conducted  under  Contract  AF  33(616)-6530  during  the  period  from 
January  1,  1963,  to  April  1,  1963,  is  reported.  This  program  is  a  research  study  to 
advance  the  state  of  the  art  of  solid  propellant  grain  design  and  consists  primarily  of 
theoretical  and  applied  analyses  of  solid  propellant  grain  designs,  considering  in  partic¬ 
ular  heat  transfer  effects,  propellant  slump,  and  the  effect  of  strain  on  burning  rate. 

An  analytical  study  is  being  undertaken  to  develop  a  two-dimensional  heat  trans¬ 
fer  analysis  for  solid  propellant  motors  subjected  to  a  nonlinear  external  thermal  environ¬ 
ment.  The  application  of  this  analysis  to  the  Advanced  Grain  Design  Computer  Program, 
previously  developed  under  this  contract,  is  also  being  investigated  to  determine  thermal 
effects  upon  ballistic  characteristics. 

The  programming  of  the  two-dimensional  heat  transfer  analysis  was  completed, 
and  several  test  cases  were  evaluated.  A  severe  stability  requirement  is  imposed  on  the 
solution  resulting  from  the  consideration  of  a  thin  case  wall  having  a  relatively  higher 
thermal  diffusivity. 

The  input  and  output  sections  of  the  Advanced  Grain  Design  Computer  Program 
were  modified,  and  certain  changes  involved  with  the  application  of  a  variable  grain  tem¬ 
perature  profile  to  the  program  were  installed. 

The  effect  of  strain  on  the  burning  rates  of  several  propellants  was  investigated. 
This  effect  increases  as  the  deviation  from  incompressible  behavior  increases. 

The  analysis  of  propellant  slump  involved  a  preliminary  analysis  of  a  finite 
length,  cylindrical  grain  having  a  spherical  head  end  and  under  acceleration  in  the  longi¬ 
tudinal  direction.  Boundary  conditions  were  considered  and  the  relevant  finite  difference 
equations  were  developed. 
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SUMMARY 


Contract  AF  33(616)-6530  is  directed  toward  the  development  of  new  and  unique 
methods  for  the  design  and  evaluation  of  solid  propellant  grains.  The  studies  undertaken 
include  theoretical  and  applied  analyses  of  solid  propellant  grain  design  from  the  stand¬ 
point  of  ballistic  and  physical  properties  of  propellants. 

The  current  program  is  divided  into  three  areas  of  study:  (1)  a  heat  transfer 
analytical  study,  (2)  the  effect  of  strain  on  burning  rate,  and  (3)  propellant  slump  analysis. 

The  major  accomplishments  during  the  period  from  January  1,  1963,  to  April  1, 
1963,  in  each  of  these  areas  are  summarized  in  the  following  three  sections. 

L _ HEAT  TRANSFER  ANALYSIS 

The  two-dimensional  heat  transfer  analysis  was  completely  programmed.  Pre¬ 
liminary  evaluation  of  the  program  indicated  that  thin  case  walls  impose  severe  stability 
requirements  upon  the  analysis  and  result  in  extremely  long  computation  times.  The  re¬ 
sults  for  two  center-perforated  grains  are  presented.  Some  modifications  to  the  program 
have  been  made  to  correct  an  error  in  the  finite  difference  equations  to  the  temperature  at 
an  internal  surface  point  on  an  irregular  boundary. 

The  Advanced  Grain  Design  Computer  Program  input  data  generator  was  modified 
to  increase  accuracy,  simplify  the  required  input  data,  and  make  the  program  more  fully 
"fail  proof."  The  output  was  revised  and  now  provides  check  results. 

The  analysis  of  solid  propellant  rocket  motors  having  a  variable  grain  tempera¬ 
ture  profile  has  been  applied  to  the  Advanced  Grain  Design  Computer  Program.  Although 
only  one  half  star  point  of  the  grain  may  be  presently  analyzed,  this  is  sufficient  for  de¬ 
bugging  purposes. 

vii 


II. 


EFFECT  OF  STRAIN  ON  BURNING  RATE 


Four  propellants  were  investigated  by  determining  the  effect  of  strain  on  Poisson's 
ratio  and  then  relating  that  parameter  to  any  changes  in  burning  rate.  There  appears  to  be 
a  correlation  between  the  deviation  from  incompressible  behavior  and  the  effect  of  strain  on 
burning  rate.  This  could  permit  a  prediction  of  the  effect  of  strain  on  burning  rate  by  an 
examination  of  Poisson's  ratio  at  that  strain. 

III.  PROPELLANT  SLUMP  ANALYSIS 

The  specification  of  the  finite  difference  equations  and  boundary  conditions  asso¬ 
ciated  with  the  elastic  stress  function  analysis  of  an  axisymmetrical  cylindrical  propellant 
charge  that  is  capped  by  a  spherical  head  end  was  accomplished.  It  is  assumed  that  the 
grain  is  under  a  constant  body  force  in  the  longitudinal  direction  and  that  the  grain  is  bonded 
to  a  rigid  case. 
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I. 


INTRODUCTION 


This  is  the  third  quarterly  report  submitted  in  partial  fulfillment  of  continuation 
studies  under  Contract  AF  33(616)-6530,  Supplement  No.  6,  Project  3059,  Task  30531, 

"A  Research  Study  to  Advance  the  State  of  the  Art  of  Solid  Propellant  Grain  Design, "  and 
summarizes  the  results  of  work  performed  during  the  period  from  January  1,  1963,  to 
April  1,  1963.  Effort  on  this  contract  supplement  was  initiated  on  July  1,  1962,  and  is 
scheduled  to  terminate  on  May  31,  1963. 

The  objective  of  this  program  is  to  advance  the  state  of  the  art  of  solid  propel¬ 
lant  grain  design  from  a  mechanical  as  well  as  ballistic  standpoint. 

Previous  studies  under  this  contract  have  included  the  development  of  a  gener¬ 
alized  computer  program  to  evaluate  the  ballistics  of  any  internal  burning,  singly  con¬ 
nected  propellant  grain  including  consideration  of  head-  and  aft-end  effects  on  the  burning 
surface  area.  The  effect  of  propellant  defects  on  ballistic  performance  was  studied  using 
windowed  test  motors  and  bombs  containing  propellant  with  artificially  prepared  voids  and 
cracks.  A  strain  analysis  was  performed  of  solid  propellant  motors  under  steady-state 
thermal  loading  employing  the  PhotoStress  technique.  Experimental  studies  were  also 
made  to  investigate  the  feasibility  of  using  an  analogous  electrical  network  in  the  solution 
of  solid  propellant  heat  conduction  problems. 

The  current  program  has  been  divided  into  three  areas  of  investigations: 

(1)  Heat  Transfer  Analysis 

(2)  Effect  of  Strain  on  Burning  Rate 


(3)  Propellant  Slump  Analysis 
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II,  HEAT  TRANSFER  ANALYSIS 

A.  Two-Dimensional  Heat  Transfer  Analysis 

An  analytical  study  is  being  performed  under  this  phase  to  consider  applications 
and  the  correlation  of  the  heat  transfer  analysis  with  the  Advanced  Grain  Design 
Computer  Program  to  determine  the  characteristics  of  solid  propellant  motors  under 
nonlinear  external  heating  conditions. 

Principal  effort  this  quarter  has  been  directed  toward  completion  of  the 
program  and  the  evaluation  of  numerous  test  cases. 

The  program  has  been  written  and  consists  of  two  sections.  Section  I  reads 
the  input  data  from  tape  and  writes  the  initial  output  on  tape.  Section  II  computes  the 
temperature  profile  at  each  time  increment  and  writes  the  profile  on  tape  as  output 
for  each  of  the  specified  print  intervals. 

Several  test  cases  were  evaluated,  two  of  which  produced  results  immediately 
usable  by  the  Rocket  Propulsion  Laboratory.  Test  oases  containing  an  irregular  in¬ 
ternal  boundary  failed  and  brought  to  light  minor  errors  in  the  basic  analysis  and 
program  logic;  these  have  been  corrected. 

A  significant  result  obtained  from  the  initial  test  case  was  the  severe 
stability  requirements  imposed  on  the  solution  by  a  thin  case  wall  having  a  relatively 
higher  thermal  diffusivity  than  the  propellant.  These  stability  requirements  result  in 
small  iteration  times  and  consequently  long  computation  times.  Since  the  validity  of 
the  case  wall  and  case-propellant  temperature  relationships  was  demonstrated 
within  the  program,  the  remaining  test  cases  were  evaluated  without  a  case  material. 


1.  Test  Case  Evaluation 


In  the  evaluation  of  the  program,  several  test  cases  were  considered. 
There  was  sufficient  variation  among  the  test  cases  so  that  each  particular  routine 
within  the  program  was  evaluated.  Each  test  case  is  described  below  and  illustrated 
in  Figures  1  through  5,  respectively. 

(1)  A  center  perforated  grain  with  a  thin  steel  case  wall.  The  center 
perforation  was  extremely  small  in  order  to  approximate  an  end 
burning  grain.  The  external  boundary  conditions  were  as  follows: 

q  =  E  cos  0  =  400  cos  0  ;  0  =  0°  -  90° 

q  a  (Trad4  -  Tos4)  =  0.  3cr(1004  -  Tos4)  2^;  0  =  90°  -180° 

The  internal  boundary  condition  was  adiabatic. 

(2)  (a)  Same  as  (1)  without  the  case  wall. 

(b)  Same  as  (2)  (a)  with  the  following  external  boundary  conditions: 

'Dfii 

q  =  E  cos  0  =  120  cos  0  ;  0  =  0°  -  90° 

q  =?0  a  <Trad4  "  T0S4)  =  0.  3o  (1004  -  Tos4)  i  0  =  90°  -  180° 

(3)  A  30°  segment  of  a  star  grain  with  no  case  wall.  The  external 
boundary  condition  was: 

,  Btu  „ 

q  ~  h0  (Tqqjjy  -  Tos)  -  8.5  (580-Tos)  £[2  >  0  —  0  -  30 

The  internal  boundary  condition  was: 

q  =  hi  (Tconv  -  T0s)  =  0.  25  (580-TOS)  dr  £^2 

(4)  A  180°  segment  of  the  same  star  grain  considered  in  (3). 

The  external  boundary  condition  was: 
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FIGURB  1.  NODAL  NETWORK  AND  BOUNDARY  CONDITIONS  FOR  TEST  CASE  (1) 
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FIGURE  2.  NODAL  NETWORK  AND  BOUNDARY  CONDITIONS  FOR  TEST  CASE  (2)  (a) 
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FIGURE  3.  NODAL  NETWORK  AND  BOUNDARY  CONDITIONS  FOR  TEST  CASE  (2>  (b) 


A4872 
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A  rg  =  0.  875  in. 

Arg  =  0. 875  in. 

A  r7  =  0. 875  in. 

Ar6  =  0. 461  in. 
Arg  =  0.461  in. 

Ar4  =  0.461  in. 

A  r3  =  0.461  in. 
A  r2  =  0. 402  in. 

A  =  0.402  in. 


Propellant  Properties 
Thermal  Conductivity  =  0. 249  Btu 

Hr-Ft 

Specific  Heat  =  0. 237  Btu 

Lb-°B 

Density  =  108.84  Lb 

Ft3 


Initial  Temperature  =  460°  R 


FIGURE  4.  NODAL  NETWORK  AND  BOUNDARY  CONDITIONS  FOR  TEST  CASE  (3) 


A4869 


-8- 


FIGURE  5.  NODAL  NETWORK  AND  BOUNDARY  CONDITIONS  FOR  TEST  CASE  (4) 
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Lifii 

q  =  E  cos  0  =  400  cos  0  ^  ^2;  0  =  0  -  90° 

q  =?0  a  (Trad4  -  Tos4)  =  0.  3  a  (1004  -  Tos4)  0  =  90°  -  180° 

The  internal  boundary  condition  was  adiabatic. 

2.  Discussion 

Test  case  1  served  to  evaluate  the  external  surface  and  case-propellant 
interface  nodal  point  relationships;  but  due  to  the  excessive  running  time  caused  by 
the  severe  stability  requirements,  no  further  information  was  obtained.  Specifically, 
the  high  thermal  diffusivity  and  small  thickness  of  the  case  wall  in  relation  to  the 
propellant  required  an  extremely  small  time  increment  to  maintain  a  stable  solution. 

Minimum  stability  is  maintained  by  using  a  time  increment  that  permits 
the  coefficient  of  the  local  nodal  point  temperature  (Tj  k)  to  be  greater  than  zero.  Thus, 
the  stability  equations  for  the  outer  surface  and  propellant  internal  nodes,  respectively, 
are  given  as  follows: 

Outer  Surface  Nodes: 

1  -  ’ 

p  Cp 


Kj,k  -  Kj-i,k 

2  A  r  o 


Kj,k  Ki,k 


Stability  restrictions  for  the  propellant  case  interface  nodes  fall  between 
the  above  two  relationships  and  should  present  no  problem  when  the  maximum  of  these 
is  used. 


For  test  case  1,  the  ratio  of  the  time  increment  required  by  the 
propellant  interior  points  to  that  required  by  the  external  surface  points  was  approximately 


In  an  attempt  to  reduce  the  computation  time  and  retain  the  case 
material,  the  following  procedure  was  outlined  and  included  within  the  program. 

A  basic  or  minimum  time  increment  is  selected,  based  upon  the  outer 
surface  stability  requirement.  Iterations  are  performed  about  the  case  for  the  outer 
surface  and  case-propellant  interface  points  alone.  This  is  continued  until  the  number 
of  iterations  about  the  case  times  the  basic  time  interval  is  less  than,  or  equal  to, 
the  time  interval  satisfying  the  internal  point  stability  requirements.  Therefore,  if 
the  basic  time  interval  is  the  time  interval  for  the  internal  point  iterations  is 


NA0g,  where  N  is  the  number  of  iterations  through  the  outer  surface  and  interface 
point. 

In  this  manner,  the  total  number  of  points  considered  is  reduced  since 
only  the  case  is  being  considered,  stability  is  satisfied  at  all  points,  and  no  loss  in 
accuracy  of  results  is  imposed.  Although  some  decrease  in  computational  time  was 
realized,  a  significant  reduction  was  not  realized  by  this  technique.  Consequently, 
further  test  cases  were  evaluated  without  a  case  wall,  which  constitutes  an  assumption 
of  a  thermally  thin  case  and/or  a  case  with  material  properties  closely  resembling 
that  of  the  propellant.  Inasmuch  as  the  case  wall-propellant  interface  calculations 
had  been  effectively  checked  out,  this  is  sufficient  for  test  purposes. 

Test  cases  (2)  (a)  and  (2)  (b)  were  successfully  evaluated,  and  the 
resulting  temperature  distributions  after  12  hours  at  the  specified  environment  are 
shown  in  Figures  6  and  7.  Intermediate  profiles  corresponding  to  every  0.  5  hours 
were  also  obtained. 

Upon  attempting  to  evaluate  test  case  3,  difficulty  was  encountered  in 
computing  the  internal  boundary  temperatures.  The  routines  within  the  program  and 
the  corresponding  finite  difference  equations  used  to  compute  these  points  were  re¬ 
viewed.  A  discrepancy  was  noted  in  the  manner  in  which  the  basic  heat  conduction 
equation  was  modified  to  account  for  an  iteration  entering  propellant  on  constant  v 
and  entering  propellant  on  constant  0  across  an  irregular  internal  boundary. 

These  modifications  were  corrected  to  the  following: 
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FIGURE  6.  TEMPERATURE  DISTRIBUTION  THROUGH  GRAIN  AFTER  12  HOURS 

FOR  TEST  CASE  2(a) 
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(1)  Entering  propellant  on  constant  0^*.  Replace  back¬ 
ward  differences  in  r  as  follows: 


Tj,k  ~  Tj-1, k  =  (T 

Ki,k 


A  rn 


(Tj,k  "  Tc) 


1  -1/2 


1  + 


rj2 


A  r 

—f 


and  Ar2  -  Kj_1>  k  -  Kj  +1>  k 


(2)  Entering  propellant  on  constant  vz£>  ,  Replace  back¬ 
ward  differences  in  0  as  follows: 


Tj,k  ~  Tj,k-1  =  (T  _  T  . 

■j.k  <Ti.k  To> 


A  0O 


1  +r  . 


-1/2 


and  Xj.k-l  '  Kj, 


k+1 


These  corrections  within  the  program  have  been  undertaken,  and  a 
satisfactory  evaluation  of  the  remaining  test  cases  should  be  available  during  the 
early  part  of  the  next  reporting  period. 

For  the  test  cases  considered,  the  heat  transfer  program  performance 
has  been  satisfactory.  The  long  computational  time  experienced  on  the  test  case 
having  a  thin  case  wall  is  characteristic  of  a  single-channel  7070.  On  either  a 
multichannel  system  or  an  internally  faster  processing  computer,  this  running  time 
would  not  be  excessively  long.  In  reality,  the  evaluation  of  the  temperature-time 


*~S  Means  "Implies" 
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history  for  the  various  imposed  environments  would  not  be  an  imposing  problem  with 
the  current  7070  computer.  However,  for  program  development  purposes,  inclusion 
of  the  case  wall  does  result  in  inefficient  use  of  computer  time. 

It  can  be  concluded  that  the  program  is  progressing  on  schedule,  and 
evaluation  of  the  remaining  test  cases  should  be  completed  during  the  first  part  of 
the  next  quarter.  At  that  time  final  documentation  of  the  program  will  be  undertaken. 

B.  Modification  of  the  Advanced  Grain  Design  Computer  Program 

The  perfection  of  the  uniform  grain  temperature  Advanced  Grain  Design 
Analysis  is  a  continuing  effort.  During  this  reporting  period,  the  input  and  output 
sections  of  the  program  were  reworked. 

The  input  data  generator  was  modified  to  increase  accuracy,  simplify  the 
required  input  data,  and  make  the  program  more  fully  "fail  proof.  "  Figure  8 
illustrates  the  generalized  rocket  motor  parameters  required  by  the  new  input  data 
generator.  The  nomenclature  for  this  program  is  listed  in  Table  I. 

A  complete  rocket  motor  test  case  was  evaluated.  The  motor  is  shown 
schematically  in  Figure  9,  while  Tables  II  through  XI  are  reproductions  of  the  output. 
These  tables  are  described  in  the  following  paragraphs. 

Table  II  -  The  first  page  of  computer  output  simply  lists  the  input  data 
exactly  as  it  appears  on  the  input  data  cards,  but  with  variable  names  printed  above 
each  word. 

Tables  III,  IV,  and  V  are  what  is  called  the  "Tape  3"  output.  Tape  3  is  the 
magnetic  tape  that  contains  the  coordinates  and  angles  between  normals  of  the  corner 
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TABLE  I 

NOMENCLATURE* 


A  -  Relative  to  Computer  Output 

I-XX  Moment  of  inertia  about  X  axis  =  /  Y2  +  Z  dm 

2  2 

I-YY  Moment  of  inertia  about  Y  axis  =  /  X  +  Z  dm 

2  2 

I-ZZ  Moment  of  inertia  about  Z  axis  =  /  X  +  Y  dm 

I-XY  "XY"  product  of  inertia  =  -  /  XYdm 

I-XZ  "XZ"  product  of  inertia  =  -  /  XZdm 

I-YZ  "YZ"  product  of  inertia  =  ~  f  YZdm 

X  A  Coordinates  of  cartesian  frame 


Y  V  oriented  (as  shown  in  Figure  10)  with 

Z  J  Z  corresponding  to  the  motor  axis 

R  Coordinates  of  polar  frame 

9  \  oriented  relative  to  XYZ  above 

Z  ,  through:  R  =-/x2  +  Y2  ;  9  =  Tan'1  (Y/X)  Z  =  Z 

S-Y  Angles  between  normals  to 

S-H  >  the  initial  burning  surface 

H-S  J  as  shown  in  Figures  10  and  11 


♦All  linear  dimensions  are  in  inches, 
All  angular  dimensions  are  in  degrees. 
All  weights  are  in  lbs. 


NOMENCLATURE 


B  -  Relative  to  Computer  Input 

(These  symbols  are  also  identified  in  Figure  8) 

Rc^  Ellipse  ratio  of  case  head  ellipse 

Zcjj  Z  translation  of  center  of  case  head  ellipse 

Acji  Major  axis  of  case  head  ellipse 

Rca  Ellipse  ratio  of  case  aft  ellipse 

Zca  Z  translation  of  center  of  case  aft  ellipse 

Aca  Major  axis  of  case  aft  ellipse 

Dport  Diameter  of  igniter  port 

Zaft  Distance  from  Z  =  0  to  aft  cutback 

Rl,  r2,  r3  Fillet  radii  in  grain  cross  section 

T-^  Least  port  diameter 

T 2  Greatest  port  diameter 


N 


Number  of  star  points 


A  6228 
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Aft  Closure: 
Elipse  or  Rati 


(1)  Constant  Distance  from  Case  Head 

(2)  Ellipse  of  Ratio  Rcore>  centered  at 
zcore  with  Major  Axis  •  Acore 

(3)  As  a  General  Conic  Section 

AR2  +  BRZ  +  CZ2  +  DR  +  EZ  +  F  «  0 


FIGURE  8.  SCHEMATIC  REPRESENTATION  OF  CROSS  SECTIONS  OF  GENERALIZED  SOLI! 

PROPELLANT  ROCKET  AMENABLE  TO  SOLUTION  WITH  AGDA  SHOWING 
REQUIRED  INPUT  PARAMETERS 


FIGURE  9.  SCHEMATIC  OF  SOLID  PROPELLANT  ROCKET  MOTOR 
FOR  WHICH  FULL  COMPUTER  OUTPUT  IS  PRESENTED  IN  TABLES  II -XI 


GENERALIZED  GRAIN  DESIGN  PROGRAM 

GRAIN  DESIGN  INPUT  GENERATOR 

COMPLETE  SOLID  PROPELLANT  ROCKET  MOTOR  EVALUATION  TEST  CA >E 

INPUT 


ZNPCT 

.0 

•  ZENO 

*  ZRND 

CONTROL  PARAMETERS 
*  ZNCORE  * 

AC 

,  .025 

*  RC 

1.4 

*  ZC 

15.7 

CASE  PARAMETERS 
*  1)1  A.M  »  AC  A 

44.05  22.025 

•  RCA 

1.4 

• 

ZC  A 

R1 

.45 

*  R2 

0.1 

CROSS  SECTION  PARAMETERS 
*  R  3  *  T 1  *  T  3  »  SN 

0.75  13.40C5  30.85  8.0 

» 

THET1C 

‘.0 

*  THET1M 
18.0 

*  THET20 
11.0 

*  THET2M  *  T  h  E  T  3  D 
37.0  2.0 

*  THt  T  3  *• 
29.5 

» 

CASEM 
i  5.76 

*  CASECG 
-0.66 

*  PROPD 
0.0632 

»  PORT  *  AETCUT 

10. J  9.733 

*  HWt  3 

6.6 

» 

DELX 
.  3 

*  DELY 

0.3 

*  DELZ 

0.5 

*  DEL  T HD  *  DEL  THM 

5.0 

• 

TABLE  REPRODUCTION  OF  INPUT  DATA  TO  PROGRAM  AGDA 
FOR  ROCKET  SHOWN  IN  FIGURE  2 
(Blank  spaces  indicate  zero) 
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GENtRAL l  Z  ED  GRAIN  DESIGN  PROGRAM 
COMPLETE  SOLID  PROPELLANT  ROCKET  MOTOR  EVALUATION  TEST  CASE 
CHECK  OUTPUT  AND  CORE  ROUND  CORRECTION 


X 

Y 

Z 

R 

THETA 

I  St  Y ) 

(lit  7> 

(  S  t  H  ) 

0.000000 

0.000000 

24.832141 

O.OOCOOO 

22.5000 

90.0000 

0.0000 

90.0000 

0. 300000 

0.124264 

24.829965 

0.324718 

22.5000 

90.0C00 

n. 7676 

90.0000 

0. 600000 

0.248528 

24.823436 

0.649435 

22.5000 

90.0000 

1.5363 

90.0000 

0.724264 

0.300000 

24.819454 

0.783938 

22.5000 

90.0000 

1.8548 

90.000.) 

0.900000 

0.372792 

24.812545 

0.974153 

22.5000 

90.0000 

2. 3059 

90.0000 

1.200000 

0.497056 

24.797279 

1.298871 

22.5000 

90.0000 

3.0776 

90.0000 

1.448528 

0.600000 

24.7813C6 

1.567876 

22.5000 

9C.OOCO 

3.7189 

90.0000 

1.500000 

0.621320 

24.777620 

1.623588 

22.5000 

90.0000 

j  .  8  5 1 9 

90.0000 

1.800000 

0.745584 

24.753543 

1.948306 

22.5000 

90.0000 

4.6294 

90.0000 

2.100000 

0.869848 

24.725020 

2.273024 

22.5000 

90.0000 

5.4110 

90.0000 

2.172792 

0.900000 

24.717425 

2.351813 

22.5000 

90.0000 

6.6013 

90. C 000 

2.400000 

0.994113 

24.692014 

2.597741 

22.5000 

90.0000 

6. 1972 

90.0000 

2.700000 

1.118377 

24.654485 

2.922459 

22.5000 

90.0000 

6.9889 

90. 0000 

2.897056 

1.200000 

24.627350 

3.135751 

22.5000 

90.0000 

7.5122 

90.0000 

3.000000 

1.242641 

24.612385 

3.247176 

22.5000 

90.0000 

7.7867 

90.0000 

3.300000 

1.366905 

24.565657 

3. 571 894 

22.5000 

90.00l:0 

8.5916 

9C.0000 

3.600000 

1.491169 

24.514243 

3.896612 

22.5000 

90.0000 

9.4042 

90.0000 

3.621320 

1.500000 

24.510409 

3.919689 

22.5000 

90.0000 

9.4623 

90.0000 

3.900000 

1.615433 

24.458071 

4.221330 

22.5000 

90.0000 

10.2256 

90.001-0 

4.200000 

1.739697 

24.397063 

4.546047 

22.5000 

90. 00r  0 

11.0566 

90.0000 

4.345585 

1.800000 

24.365692 

4.703627 

22.5000 

99.0COO 

11.4637 

90.0000 

4.500000 

1.863961 

24.331139 

4.870765 

22.5000 

90.0000 

11.8932 

90.0000 

4.800000 

1.988225 

24.260193 

5.195482 

22.5000 

9C.000 0 

12.7516 

90.0000 

4.069849 

2.100000 

24.192010 

5.487565 

22.5000 

90.0000 

13.5301 

90.  1.1COO 

4. 100000 

2. 112489 

24. 1841 30 

5.520200 

22.5000 

90.0000 

l  3.6173 

9 J. OOOd 

5.400000 

2.236753 

24. 102825 

5.844918 

22.5000 

90.0000 

14.4981 

90.1,000 

4.700000 

2.361017 

24.016152 

6. 169636 

22.5000 

90.0000 

15. 3937 

90.0000 

5.794113 

2.400000 

23.987831 

6.271502 

22.5000 

90.0000 

15.6781 

90.  '.000 

6.000000 

2.485281 

23.923967 

6.494353 

22.5000 

90.0005 

16.306) 

90.0000 

6.190224 

2.564077 

23.862589 

6. 700251 

22.5000 

67.5003 

16.8941 

73.1059 

6.200973 

2.540124 

23.862342 

6.701067 

22.2756 

64.157C 

16.8965 

73.1373 

4.213159 

2.516738 

23.861594 

6.703530 

22.0512 

6C.7989 

16.9035 

73.2318 

6.226756 

2.494008 

23.860341 

6.707649 

21.8276 

57.4261 

16.9154 

73.3899 

6.241735 

2.472022 

23.858581 

6.713430 

21.6059 

54.0383 

16.9320 

73.6118 

6.258063 

2.450867 

23.856312 

6. 720871 

21.3870 

50.6354 

16.9535 

73.8974 

6.275703 

2.430629 

23.853541 

6. 729963 

21.1718 

47.2166 

13.979a 

74.246  3 

6.294610 

2.411390 

23.850258 

6. 740690 

20.9612 

43.7816 

1 7.0103 

74.6  >80 

6.306948 

2.400000 

23.847974 

6.748154 

20.8335 

41.6435 

17.0321 

74. 9444 

6.313560 

2.394236 

23.846708 

6.752288 

20.7678 

40.5266 

17.0440 

7  >.  103u 

6.333531 

2.378125 

23.842719 

6. 765286 

20.5802 

37.2590 

17.0313 

75.601 7 

6.354496 

2.3631C6 

23.838291 

6.779667 

20.3991 

33.9730 

17. 1 2  3  f. 

76. 1636 

TABLE  Hi  REPRODUCTION  OF  "TAPE  3"  AGDA  OUTPUT 
(Figures  3  and  4  illustrate  meaning  of  these  numbers) 
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GENERALIZED  GRAIN  DESIGN  PROGRAM 


COMPLETE  SOLID  PROPELLANT  ROCKET  MOTOR  EVALUATION  TEST  CASE 


X 

Y 

Z 

R 

THETA 

(S,Y) 

«H,Z) 

1  S»  H) 

6. 376430 

2.349228 

23.833432 

6.795420 

20.2250 

30.6698 

17.1686 

76.7583 

6.399276 

2.336561 

23.828146 

6.812507 

20.0586 

27.3432 

17.2180 

77.4144 

6.422972 

2.325173 

23.822445 

6.830885 

19.9C07 

23.9954 

17.2712 

78. 1206 

6.447448 

2.315128 

23.816332 

6.850504 

19.7520 

20.6262 

17.3281 

78.8749 

6.472635 

2.306491 

23.809828 

6.871310 

19.6133 

17.2355 

17. 3986 

79.6755 

6.498455 

2.299316 

23.802949 

6.893241 

19.4851 

13.8230 

17.4523 

80.5205 

6.524826 

2.293660 

23.795703 

6.916230 

19.3680 

10. 3884 

17.5194 

81.4076 

6.551664 

2.289572 

23.788119 

6.940205 

19.2627 

6.9313 

17.5893 

82.3344 

6.578879 

2.287099 

23.7802 14 

6.965089 

19.1696 

3.4515 

17.6620 

83.2934 

6.606375 

2.286284 

23.772009 

6.990800 

19.0893 

-0.0512 

17.7373 

84.2968 

6.623244 

2.286614 

23.766868 

7.006851 

19.0468 

-2.1993 

17.7844 

84.9216 

6.923244 

2.298138 

23.672172 

7.294707 

18.3633 

-2. 1998 

13.6373 

84.3961 

7.223244 

2.309663 

23.572312 

7.583522 

17.7319 

-2.1998 

19.5112 

34.8689 

7.523244 

2.321188 

23.467124 

7.873190 

17. 1469 

-2.1998 

20.4075 

34. 8401 

7.823244 

2.332713 

23.356410 

8.163620 

16.6034 

-2.1998 

21.3273 

87.8093 

8.123244 

2.344238 

23.239977 

8.454735 

16.0973 

-2.1993 

22.2739 

84.7757 

8.423244 

2.355763 

23.117592 

8.746465 

15.6249 

-2.1993 

23.2479 

87.7420 

8.723244 

2.367287 

22.988992 

9.038752 

15.1831 

-2.1998 

24.2522 

64.7051 

8. 790053 

2.369846 

22.959484 

9.103912 

15.0885 

-2.2002 

24.4802 

34.6968 

8.796120 

2.369896 

22.956809 

9.109782 

15.0789 

1.2763 

24.5008 

83.2938 

8.802218 

2.369575 

22.954161 

9.115587 

15.0670 

4.7762 

24.5212 

81.9012 

3. 810253 

2.368580 

22.950736 

9.123088 

15.0478 

9.4166 

24.5476 

80.0930 

9.110253 

2.318826 

22.821029 

9.400727 

14.2802 

9.4166 

25.5395 

80.0219 

9.410253 

2.269072 

22.684520 

9.679956 

13.5568 

9.4166 

26.5690 

79.9461 

9.710253 

2.219317 

22.540878 

9.960642 

12.8741 

•5.4166 

27.6394 

79.8653 

10.010253 

2.169563 

22.389695 

10.242663 

12.2288 

9.4166 

28.7544 

79. 7791 

10.310253 

2.119809 

22.230522 

10.525915 

11.6182 

9.4166 

29.9135 

79.6870 

10.429692 

2.100000 

22.164814 

10.639007 

11.3842 

9.4166 

30.3967 

79.6485 

10.610253 

2.070054 

22.062839 

10.810299 

11.0397 

9.4166 

31.1369 

79.5884 

10.910253 

2.020300 

21.886056 

11.095730 

10.4909 

9.4166 

32.4158 

79.4323 

11.210253 

1.970546 

21.699459 

11.382127 

9.9696 

9.4166 

33.7628 

79.3695 

11.510253 

1.920792 

21.502233 

11.669419 

9.4740 

9.4166 

35. 1866 

79.2476 

11.810253 

1.871037 

21.293379 

11.957543 

9.0023 

9.4166 

36.6984 

79.1161 

12.110253 

1.821283 

21.071704 

12.246439 

8.5527 

9.4166 

38.3119 

78.9738 

12.238580 

1. 800000 

20.972610 

12.370240 

8.3668 

9.4166 

39.0372 

78. 9093 

12.410253 

1.771529 

20.835731 

12.536055 

8.1239 

9.4166 

40.0441 

78.  8194 

12.710253 

1.721774 

20.583648 

12.826341 

7.7145 

9.4166 

41.9170 

78.6511 

13.010253 

1.672020 

20.313044 

13.117253 

7.3233 

9.4166 

43.9598 

70.4663 

13.310253 

1.622266 

20.020780 

13.408749 

6.9490 

9.4166 

46.2119 

7R.2637 

13.610253 

1.572511 

19.702569 

13.700794 

6.5907 

9.4166 

43.7283 

78. C  332 

13.910253 

1.522757 

19.352159 

13.993352 

6.2473 

9.4166 

51.5909 

77.7858 

TABLE 


REPRODUCTION  OF  PAGE  2  OF  "TAPE  3"  AG  DA  OUTPUT 
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GENERAL12ED  GRAIN  DESIGN  PROGRAM 
COMPLETE  SOLID  PROPELLANT  ROCKET  MOTOR  EVALUATION  TEST  CASE 


X 

Y 

Z 

R 

THETA 

t  S»  Y) 

<H,2) 

(S,b) 

14.047468 

1.500000 

19.178724 

14.127326 

6.0950 

9.4166 

53.0480 

77.6597 

14.210253 

1.473003 

18.959952 

14.286392 

5.9180 

9.4166 

54.9230 

77.5000 

14.510253 

1.423248 

18.509560 

14.579885 

5.6020 

9.4166 

58.9621 

7  7.1710 

14.783834 

1.377875 

18.021C92 

14. 347904 

5.3247 

9.4166 

63.6179 

74. J234 

14. 828248 

1.369126 

17.931825 

14.891320 

5.2753 

12.0753 

64.5032 

73.6696 

14.872167 

1.357677 

17.8404*72 

14.934009 

5.2161 

16.3431 

65.4207 

70.4782 

14.915425 

1.343554 

17.747043 

14.975814 

5.1472 

19.8200 

66.3712 

67.2500 

14.957859 

1.326786 

17.651610 

15.016587 

5.0690 

23.3062 

67.  3550 

63.9351 

14.999306 

1.307408 

17.554289 

15.056177 

4.9016 

26.8020 

63.3716 

60.603d 

15.039606 

1.285478 

17.455135 

15.094442 

4.8854 

30.3076 

69.4218 

57.3465 

1 5.078603 

1.261048 

17.354268 

15.  131242 

4.7806 

33.8236 

70.5046 

53.9730 

15.116142 

1.234188 

17. 2517C2 

15.166441 

4.6677 

37.3504 

71.6207 

5  0 . 3  j  3  6 

15.152073 

1.204966 

17.147496 

15.199909 

4.5469 

40.8891 

72.7701 

47.1178 

15.157753 

1.200000 

17.130406 

15.205179 

4.5265 

41.4655 

72.9600 

46. 5542 

15.186250 

1.173473 

17.041721 

15.231520 

4.4186 

44.4400 

73.9524 

43.6357 

15.218533 

1.139780 

16.934473 

15.261154 

4.2831 

48.0054 

75.1668 

40.1153 

15.247043 

1.106174 

16.832222 

15.287116 

4.1495 

51.3725 

76.339" 

36.7720 

15.273603 

1.070854 

16.728827 

15.311095 

4.0105 

54.7490 

77.5379 

33.3934 

15.298086 

1.033945 

16.624559 

15. 332986 

3.8665 

58.1331 

78.7600 

30.  O'*  18 

15.320389 

0.995584 

16.519368 

15.352703 

3. 7181 

61.5234 

30.0054 

26.5813 

15.340396 

0.955933 

16.413635 

15.370150 

3.5658 

64.9168 

31.2683 

23. 1432 

15.357996 

0.915242 

16.307463 

15. 385242 

3.4104 

63.3041 

82.5481 

19. 697  4 

15.363869 

0.900000 

16.268152 

15.390206 

3.3525 

67.5520 

83.024? 

13.4245 

15.373234 

0.873366 

16.200389 

15.398021 

3.2515 

71.7069 

83.8476 

16.2224 

15.386031 

0.830378 

16.092272 

15.408421 

3.0892 

75.1359 

33.1679 

12.7148 

15.396291 

0.786414 

15.982956 

15.416361 

2.9240 

78.5863 

86.5094 

9.1753 

15.403937 

0.741617 

15.872899 

15.421778 

2.7564 

82.0575 

87.3650 

5.607? 

15.408868 

0.696139 

15.762443 

15.424584 

2.5867 

85.5527 

89.2285 

2.0140 

15.410416 

0.670597 

15.707111 

15.424999 

2.4917 

87.5083 

39.9121 

0.0369 

15.413326 

0.600000 

15.707111 

15.424999 

2.2292 

87.7707 

89.9121 

0 . 08o  9 

15.422082 

0.300000 

15.707111 

15.424999 

1.1144 

88.8856 

89.9121 

3.0869 

15.423409 

0.221540 

15.707111 

15.424999 

0.8229 

89.1771 

39.9121 

0.0850 

15.425000 

0.000000 

15.707111 

15. 424999 

0.0000 

90.0000 

89.9121 

C.C383 

TABLE  REPRODUCTION  OF  PAGE  3  OF  AGDA  COMPUTER  OUTPUT 


GENERALIZED  GRAIN  DESIGN  PROGRAM 
COMPLETE  SOLID  PROPELLANT  ROCKET  MOTOR  EVALUATION  TEST  CASE 

SUPPLEMENTARY  CHECK  RESULTS 


TOTAL 

HEAD  END 

TOP  SECT. 
ST.  THRU 

CENTER  SECT. 

ST.  THRU 

AFT  END 

SURFACE 

4022.39 

509.56 

878.05 

2634.79 

o 

• 

o 

o 

VOLUME 

9434.28 

2899.71 

6534.57 

0.00 

PORT  AREA 

*  416.0262 

TABLE  VI.  CHECK  RESULTS 
(These  results  are  based  only  upon  "Tape  3"  data  and 
therefore  serve  to  validate  the  majority  of  AGDA. ) 


COMPUTEO  TRACE  OF  INTERNAL  GEOMETRY 


COMPLETE  SOLID  PROPELLANT  ROCKET  MOTOR  EVALUATION  TEST  CASE 


WEB 

CHAMBER 

CHAMBER 

FIRST  l 

BURNOUT 

SURFACE 

VOLUME 

MGMENT 

(PERCENT) 

(SQ-IN) 

( CU-IN) 

( IN**4 ) 

0 

5.5747E 

03 

1.7151E 

04 

3.8688b 

05 

5 

5.7019E 

03 

1.9012E 

04 

3.7302E 

05 

10 

5.6624E 

03 

2.0901E 

04 

3.5862b 

05 

15 

5.5267E 

03 

2.2747E 

04 

3.44206 

05 

20 

5.4046E 

03 

2.4546E 

04 

3.2976b 

05 

25 

5.4202E 

03 

2.6331E 

04 

3.1509b 

05 

30 

5.4407E 

03 

2.8123E 

04 

3. 000 lr 

05 

35 

5.4533E 

03 

2.9920E 

04 

2.84496 

05 

40 

5.4620E 

03 

3.1720E 

04 

2.6350ir 

05 

45 

5.4666E 

03 

3.3522E 

04 

2.5201b 

05 

50 

5.4663E 

03 

3.5325E 

04 

2.34996 

05 

55 

5.4621E 

C3 

3.7127E 

04 

2.1741b 

05 

60 

5.45246 

03 

3.8926E 

04 

1 .992  IE 

05 

65 

5.4373E 

03 

4.0722E 

04 

1.8039E 

05 

70 

5.4150E 

03 

4.2511E 

04 

1.6089b 

05 

75 

5.3858E 

03 

4.4292E 

04 

1.4066E 

05 

80 

5.3518E 

03 

4.6060E 

04 

1.19716 

05 

85 

5.3232E 

03 

4.7818E 

04 

9.7932E 

04 

90 

5.30156 

03 

4.9568E 

04 

7.53056 

04 

95 

5.2695E 

03 

5.1309E 

04 

5.17996 

04 

100- 

5.2273E 

03 

100+ 

4.4455E 

03 

5.3027E 

04 

2.75636 

04 

105 

2.4510E 

03 

5.4002E 

04 

1.73046 

04 

110 

1.8014E 

03 

5.4692E 

04 

1.0595E 

04 

115 

1.2943E 

03 

5.5193E 

04 

5.9546E 

03 

120 

8.7868E 

02 

5.5546E 

04 

2.83866 

03 

125 

5.0963E 

02 

5.5771E 

04 

9.6954E 

02 

130 

1.6718E 

02 

5.5878E 

04 

1.2255b 

02 

WEB  THICKNESS  = 

6. 

600 

HEIGHT  = 

15. 

700 

DIAMETER  = 

44. 

050 

CASE  VOLUME  = 

55894. 

089 

INHIBITED  PORT  DIAMETER  = 

1C. 

000 

TABLE  VII.  FIRST  PAGE  OF  ORIGINAL  PROGRAM  AG  DA  OUTPUT 


COMPUTED  TRACE  OP  MOMENTS  OF  INERTIA 


COMPLETE  SOLID  PROPELLANT  ROCKET  MOTOR  EVALUATION  TEST  CASE 


WEB 

BURNOUT 

(PERCENT) 

I-XX 

( IN**5) 

I-YV 

(  IN»*5) 

1-2/ 

< I N*«5 ) 

0 

1.3464E 

07 

1.34646 

07 

1.05346 

07 

5 

1.3054E 

07 

1.30546 

07 

1.02806 

07 

10 

1.2619E 

07 

1.2619E 

07 

1 . 00006 

07 

15 

1.2171E 

07 

1.21716 

07 

9.70076 

06 

20 

1.1710E 

07 

1.17106 

07 

9.38116 

06 

25 

1.1230E 

07 

1.12306 

07 

9.03946 

06 

30 

1.0726E 

07 

1.07266 

07 

8.67346 

06 

35 

1.0199E 

07 

1.0199E 

07 

8.2833.: 

06 

40 

9.6484F 

06 

9.6484E 

06 

7.869H- 

06 

45 

9.07326 

Ob 

9.0732E 

06 

7.43136 

06 

50 

8.4736E 

06 

8.4736E 

06 

6.97C16 

06 

55 

7.8491E 

06 

7.8491E 

06 

6.48596 

06 

60 

7.1994E 

06 

7.1994E 

06 

5.9794*- 

06 

65 

6.5245E 

06 

6.52456 

06 

5  •  4  5 1 3 1 

06 

70 

5.8240E 

06 

5.8240E 

06 

4.90266 

06 

75 

5.0975E 

06 

5.0975E 

06 

4.33456 

06 

80 

4.3465E 

06 

4.34656 

06 

3.74916 

06 

85 

3.5690E 

06 

3.5690E 

06 

3.14646 

06 

90 

2.7642E 

06 

2.7642E 

06 

2.5266c 

06 

95 

1.9318E 

06 

1.93186 

06 

1.8905c 

06 

100 

1 .07896 

06 

1.07896 

06 

1.24476 

06 

105 

6.89746 

05 

6.89746 

05 

8.3811c 

0  5 

110 

4.2827E 

05 

4.2827E 

05 

5.4197E 

05 

115 

2.44066 

05 

2.44066 

05 

3.21246 

05 

120 

1. 18196 

05 

1.1819E 

05 

1.62096 

05 

125 

4.0706E 

04 

4.07066 

04 

5.8056c 

04 

130 

4.04896 

02 

6.3055E 

01 

8. 5904c 

01 

TABLE  VIII.  PAGE  2  OF  ORIGINAL  PROGRAM  AG  DA  OUTPUT 


COMPUTED  TRACE  OF  PRODUCTS  OF  INERTIA 


COMPLETE 

SOLID  PROPELLANT 

ROCKET  MOTOR 

EVALUATION  TEST 

WEB 

BURNOUT 

I-XY 

I-XZ 

I-YZ 

(PERCENT) 

( IN**5) 

CIN**5) 

( IN**5 ) 

0 

8.2441E-02 

3.A223E-02 

1. 2320E-01 

5 

8.0341E-02 

3.3352E-02 

1.20076-01 

10 

7.8064E-02 

3.2A1AE-02 

1. 16696-01 

15 

7. 5633E-02 

3.  14326-02 

1.1 3156-01 

20 

7.3051E-02 

3.0A08E-02 

1 • 0947F-01 

25 

7.0301E-02 

2. 93306-02 

1.05596-01 

30 

6.7370E-02 

2.8191E-02 

1.01496-01 

35 

6.4259E-02 

2.69866-02 

9.71506-02 

AO 

6.0969E-02 

2.5711E-02 

9.25596-02 

45 

5.7502E-02 

2. A362E-02 

8.77036-02 

50 

5.3861E-02 

2. 2935E-02 

8.2567c-02 

55 

5.0046E-02 

2. 1A27E-02 

7.71376-02 

60 

4.6062E-02 

1.98326-02 

7.13966-02 

65 

A.  1912E-02 

1.81A7E-02 

6.  5329l:-02 

70 

3.7603E-02 

1.6366E— 02 

5. 8917c-02 

75 

3.31A0E-02 

1.AA83E-02 

5.21386-02 

80 

2.85386-02 

1.2A9A6-02 

4. 49796-02 

85 

2.3797E-02 

1. 0390E-02 

3.74036-02 

90 

1.8918E-02 

8.1636E-03 

2. 9389E-02 

95 

1. 3912E-C2 

5.8095E-03 

2. 09146-02 

100 

8.8328E-03 

3. 3369E-03 

1.20136-02 

105 

5. 8069E-03 

2. 12A1E-03 

7.64686-03 

110 

3.6828E-C3 

1.  3135E-03 

4.72876-03 

115 

2.1A36E-03 

7.  AA26E-0A 

2.67936-03 

120 

1.0637E-03 

3. 5736E-0A 

1.2865E-03 

125 

3.7522E-0A 

1.2278E-0A 

4.42016-04 

130 

O.OOOOE  OC 

1.5602E  02 

O.OOOOt  00 

TABLE  IX.  PAGE  3  OF  ORIGINAL  PROGRAM  AG  DA  OUTPUT 


THIOKOL-ELKTON 

ADVANCEC  GRAIN  DESIGN  COMPUTER  PROGRAM 

COMPLETE  SOL  10  PROPELLANT  ROCKET  MOTOR  EVALUATION  TEST  CASE 

THEORETICAL  ANALYSIS  BY  F.E. MOORE 
PROGRAMMING  ANALYSIS  BY  D.H.FREOERICK 


V  1.  0 

BURNING 

CHAMBER 

PROP.  PROP. 

PROPELLANT 

INERTIA 

URN 

SURFACE 

VOLUME 

VOLUME 

CG 

I-XX 

I-YY 

"C.T 

(SQ-INI 

ICU-IN) 

(CU-IN) 

(IN) 

- ( SO- 

ET«CU-I\ 

0 

5574.66 

17151.03 

38743.06 

9.986 

93500.6 

93500.6 

5 

5701.89 

19011.64 

36882.45 

10.114 

90655.6 

90655.6 

to 

5662.43 

20901.43 

34992.66 

10.248 

87634.9 

37634.5 

15 

5526.67 

22747.39 

33146.70 

10.3e4 

84521.2 

84521.2 

20 

5404.55 

24546.06 

31348.03 

10.519 

31318.8 

81318.8 

25 

5420.16 

26330.97 

29563.12 

10.658 

77984.9 

77984.9 

30 

5440.68 

28122.73 

27771.36 

10.803 

74488.9 

74488.9 

35 

5453.31 

29919.56 

25974.53 

10.953 

70828.9 

70828.9 

4  0 

5461.95 

31719.94 

24174.15 

11.107 

67002.6 

67002.6 

45 

5466.61 

33522.22 

22371.87 

11.265 

63008.6 

63008.6 

■>0 

5466.34 

35324.93 

20569.16 

11.425 

58844.7 

58844.7 

^5 

5462.09 

37126.86 

18767.23 

11.584 

54508.0 

54508.0 

oO 

5452.42 

38926.37 

16967.72 

11.741 

49996. 1 

49996.1 

65 

5437.26 

40721.74 

15172.35 

11.889 

45308.7 

45308.7 

70 

5415.02 

42510.87 

13383.22 

12.021 

40444.2 

40444.2 

73 

5385.84 

44291.97 

11602.12 

12.124 

35399.6 

35399.6 

20 

5351.81 

46059.71 

9834.38 

12.173 

30184.3 

30134. 3 

-.5 

5323.25 

47817.55 

8076.54 

12.126 

24784.6 

24784.6 

1 0 

5301.46 

99567.61 

6326.48 

11.903 

19195.9 

19195.9 

>5 

5269.46 

51309.21 

4584.88 

11.298 

13415.4 

13415.4 

1  00 

5227.27 

53026.60 

2867.49 

9.612 

7492.2 

7492.2 

'  AX 

5701.89 

53026.60 

38743.06 

12.173 

93500.6 

93500.6 

■  1 N 

5227.27 

17151.03 

2867.49 

9.612 

7492.2 

7492.2 

r.  VG 

5438.34 

35269.35 

20624.73 

11.140 

55739.0 

55739.0 

loo* 

1  >5 

4445.55 

2450.99 

54001.86 

1892.23 

4789.8 

4789.8 

1  10 

1801.36 

54691.78 

1202.31 

2974.1 

2974.1 

!  13 

1294.28 

55192.89 

701.20 

1694.9 

1694.9 

!  20 

878.68 

55545.79 

348.30 

820.8 

820.3 

125 

509.63 

55771.24 

122.85 

282.7 

282.7 

1  30 

167.18 

55878.48 

15.61 

2.8 

0.4 

PER 

CENT  PROPELLANT  (  0- 

100  PCT  WEB) 

BY  VOLUME 

=  64.18 

PER 

CENT  PROPELLANT  ( 100- 

130  PCT  WEB) 

BY  VOLUME 

=  5.13 

TABLE  A.  PAGE  1  OF  NEW  PROGRAM  AG  DA  OUTPUT 


TENSOR 

1-17 


73154.3 

71385.5 

69447.8 

67365.9 

65146.9 
6277?. 9 

60232.1 

57522.7 

54646.6 
51606.0 

42403.2 

45040.8 

41523.6 

37856.2 
34046.0 

30100.9 

26035.4 

21850.3 

17543.6 

13128.6 

8644.1 

73154.8 

8644.1 

45593.2 


5820.2 
376  5.  7 

2230.8 

1125.6 

403.2 


(The  max. ,  min. ,  and  average  values  refer  only  to  0-9100  pet  web 
The  100  +  surface  area  is  the  surface  just  after  web  burnout. ) 
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THIOKOL-ELKTON 

ADVANCED  GRAIN  DESIGN  COMPUTER  PROGRAM 

COMPLETE  SOLID  PROPELLANT  ROCKET  MOTOR  EVALUATION  TEST  CASE 

THEORETICAL  ANALYSIS  BY  F.E. MOORE 
PROGRAMMING  ANALYSIS  BY  D.H. FREDERICK 


KtB 

BURNING 

PROP. 

SYSTEM 

SYSTEM 

PROPELLANT 

INERt'I  A 

TENSOR 

i  URN 

SURFACE 

HEIGHT 

HEIGHT 

CG 

I-XX 

I  —  YY 

I -U. 

(SQ-IN) 

(LBS) 

(LBS) 

(IN) 

*  1 *Luw 1 

0 

5574.66 

2448.56 

3154.32 

7.6C4 

5909.2 

5909.2 

4623.4 

5 

5701.89 

2330.97 

3036.73 

7.610 

5729.4 

5729.4 

4511.6 

10 

5662.43 

2211.54 

2917.30 

7.609 

5538.5 

5538.5 

4389. 1 

15 

5526.67 

2094.87 

2800.63 

7.6Q1 

5341.7 

5341.7 

4257.5 

20 

5404.55 

1981.20 

2686.96 

7.583 

5139.3 

5139.3 

4117.3 

25 

5420.16 

1868.39 

2574.15 

7.555 

4928.6 

4928.6 

3967.3 

30 

5440.66 

1755.15 

2460.91 

7.515 

4707.7 

4707.7 

3806.7 

35 

5453.31 

1641.59 

2347.35 

7.461 

4476.4 

4476.4 

3635.4 

40 

5461.95 

1527.81 

2233.57 

7.389 

4234.6 

4234.6 

3453.7 

45 

5466.61 

1413.90 

2119.66 

7.294 

3982.1 

3982.1 

3261.5 

30 

5466.34 

1299.97 

2005.73 

7.172 

3719.0 

3719.0 

3059.1 

}5 

5462.09 

1186.09 

1891.85 

7.017 

3444.9 

3444.9 

2846.6 

bD 

5452.42 

1072.36 

1778.12 

6.819 

3159.8 

3159.8 

2624.3 

65 

5437.26 

958.89 

1664.65 

6.569 

2863.5 

2863.5 

2392.5 

70 

5415.02 

845.82 

1551.58 

6.253 

2556.1 

2556.1 

2151.7 

75 

5385.84 

733.25 

1439.01 

5.854 

2237.3 

2237.3 

1902.4 

MO 

5351.81 

621.53 

1327.29 

5.349 

1907.6 

1907.6 

1645.4 

35 

5323.25 

510.44 

1216.20 

4.706 

1566.4 

1566.4 

1380.9 

iO 

5301.46 

399.83 

1105.59 

3.883 

1213.2 

1213.2 

1108.9 

35 

5269.46 

289.76 

995.52 

2.820 

847.9 

847.9 

829.7 

1  00 

5227.27 

181.23 

886.99 

1.439 

473.5 

473.5 

546.3 

VAX 

5701.89 

2448.56 

3154.32 

7.610 

5909.2 

5909.2 

4623.4 

MIN 

5227.27 

181.23 

886.99 

1.439 

473.5 

473.5 

546.3 

AVG 

5438.34 

1303.48 

2009.24 

6.338 

3522.7 

3522.7 

2881.5 

i  00* 

4445.55 

105 

2450.99 

119.59 

825.35 

0.761 

302.7 

302.7 

367.8 

1 10 

1801.36 

75.99 

781.75 

0.261 

188.0 

188.0 

237.9 

115 

1294.28 

44.32 

750.08 

-0.119 

107.1 

107.1 

141.0 

120 

878.68 

22.01 

727.77 

-0.394 

51.9 

51.  9 

71.1 

125 

509.63 

7.76 

713.52 

-0.567 

17.9 

17.9 

25.5 

130 

167.18 

0.99 

706.75 

-0.648 

0.2 

0.0 

3.0 

TABLE  XL  PAGE  2  OF  NEW  PROGRAM  AGDA  FINAL  OUTPUT 

(This  page  includes  the  weight  of  inert  parts  as  well  as  propellant 
density  to  obtain  actual  weight,  CG  travel,  etc. ) 
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of  the  rocket  motor  core .  These  coordinates  are  defined  by  the  intersection  of  the  side  of 
a  star  point  and  the  core  head.  Figures  10  and  11  illustrate  the  area  just  described. 

At  the  conclusion  of  the  Tape  3  output,  the  supplementary  check  results  are 
printed  (Table  VI).  These  results  are  obtained  from  Tape  3  data  only,  and  no  use  is 
made  of  the  bulk  of  the  Advanced  Grain  Design  Analysis  program.  Initial  surface 
area  and  initial  chamber  volume  are  calculated. 

The  next  three  pages  of  output  (Tables  VII,  VIII,  and  IX)  are  in  the  original 
output  format.  Table  VI shows  burning  surface,  chamber  volume,  and  the  first  moment 
of  the  propellant  along  the  Z  axis.  Table  VIII  shows  the  moments  of  inertia, 

while  Table  IX  shows  the  products  of  inertia. 

The  final  two  pages  of  computer  output  (Tables  X  and  XI)  represent  the  new 
output  data  format.  These  pages  are  identical  except  for  the  fact  that  Table  X  is 
independent  of  propellant  density,  while  Table  XI  utilizes  this  quantity  for 
determination  of  propellant  weight,  etc.  Case  weight  and  payload  weight,  if  specified 
to  the  program,  are  used  in  conjunction  with  the  computed  propellant  center  of  gravity 
travel  to  compute  the  C.  G*  travel  of  the  entire  system. 

The  three  ^dimensional  AGDA  computer  program  appears  to  operate  properly 
for  ^11  test  cases  evaluated.  Reservatiqn  is  implied  here  since  the  rigorous  analysis  of  the 
AGDA  is  unmatched  by  manual  analysis  and  other  known  computer  analyses .  Accordingly, 
minor  differences  between  predicted  and  AGDA  determined  answers  cannot  be  readily 
resolved.  Despite  its  complexity,  the  program  is  extremely  easy  to  use.  The  input 
data  generator  has  all  but  eliminated  errors  caused  by  incorrect  input  data.  In  addition 


1 
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The  core  or 
central  perfora¬ 
tion  of  a  solid 
propellant  rocket 
motor. 


2 


Magnified  view  of 
shaded  half  star  point 
of  1.  The  dots  indicate 
points  at  which  cartesian 
and  polar  coordinates  are 
computed  for  the  Tape  3 
output. 


u  if  j  fr 


X 
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Longitudinal  cross  section  of  rocket  motor  which  reveals  the  normal  if  and 
the  angle  Z-H 


Axial  cross  section  of  one  half  star  point  which  reveals  the  normal^ and  the 
angle  S-Y 


The  angle  S-H  is  simply  the  angle  between  normals  if  and  it  Reference  back  to 
Figure  10  will  show  a  three  dimensional  illustration  of  both  if  and  if  and  the  angle 
S-H.  , 


FIGURE  11.  DEFINITION  OF-  NORMALS  f  AND  if 
TOGETHER  WITH  THEIR  ASSOCIATED  ANGLES. 
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to  an  easily  readable  format,  the  AGDA  output  provides  check  results  that  may  be 
compared  to  the  main  program  output  for  the  purpose  of  checking  program  accuracy. 

C,  Irregular  Temperature  Profile 

The  straight-through  portion  of  the  AGDA  computer  program  was  modified 
according  to  the  theory  presented  in  the  previous  quarterly  report1  in  order  to  analyze 
solid  propellant  rocket  motors  having  a  variable  grain  temperature  profile.  At  the 
present  time,  these  modifications  permit  temperature  variation  over  only  one  half 
star  point  of  the  rocket  motor.  Test  results  for  debugging  purposes  may  be  obtained 
in  a  much  shorter  period  of  time  when  only  one  half  star  point  is  analyzed.  Although 
in  principle  an  entire  rocket  motor  can  be  analyzed  by  manually  executing  one  half 
star  point  at  a  time,  the  capability  of  automatically  performing  the  complete  analysis 
will  be  included  once  the  program  is  fully  checked  out  for  the  half  star  point  case. 

Three  test  cases  have  been  run  with  the  present  version  of  the  variable 
temperature  program.  Each  test  case  involved  the  straight-through  section  of  the 
rocket  motor  shown  in  Figure  9. 

Case  I  -  Constant  temperature  —  Variable  temperature  option  disabled 

Case  II  -  Constant  temperature  —  Variable  temperature  option  used  with 

table  of  equal  temperature 

Case  III  -  Variable  temperature  option  used  with  nonuniform  grain  temperature. 

Inspection  of  Tables  XII  and  XIII,  which  are  reproductions  of  the  final  page 
of  output  from  test  case  I  and  test  case  II,  respectively,  shows  that  these  outputs  are 


identical.  This  is  the  desired  result. 
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ADVANCED  GRAIN  DESIGN  COMPUTER  PROGRAM 

STRAIGHT  THRU  SECTION  BROKEN  BACK  STAR  VARIABLE  TEMP  OPTION  DISABLED 

THEORETICAL  ANALYSIS  BY  F.E. MOORE 
PROGRAMMING  ANALYSIS  BY  D.H. FREDERICK 


WEB 

BURNING 

CHAMBER 

PROP.  PROP. 

PROPELLANT 

INERTIA 

TENSOR 

BURN 

SURFACE 

VOLUME 

VOLUME 

CG 

I-XX 

I-YY 

1-22 

PCT 

ISO-INI 

tCU-IN) 

(CU-1N) 

(IN) 

- ( SQ-FT*CU-IN) 

C 

2634.05 

6556.54 

1737C.08 

7.B50 

28072.6 

28072.6 

36323.3 

5 

2666.59 

7430.92 

16495.70 

7.850 

27138.3 

27138.3 

35452.4 

10 

2620.76 

8310.01 

15616.61 

7.850 

26163.4 

26163.4 

34505.8 

15 

2530.06 

9159.67 

14766.95 

7.850 

25173.5 

25173.5 

33495.6 

20 

2445.54 

9978.20 

13948.42 

7.850 

24171.6 

24171.6 

32425.8 

25 

2426.25 

10781.46 

13145.16 

7.850 

23145.0 

23145.1 

31289.4 

30 

2409.74 

11579.18 

12347.44 

7.850 

22085.2 

22085.2 

30080.0 

35 

2393.22 

12371.46 

11555.16 

7.850 

20991.5 

20991^5 

28796.8 

40 

2376.71 

13158.30 

10768.33 

7.850 

19863.8 

19863.8 

27439.3 

45 

2360.20 

13939.68 

9986.94 

7.850 

18701.9 

18701.9 

26007.2 

50 

2343.68 

14715.61 

9211.01 

7.850 

17505.7 

17505.7 

24500.3 

55 

2327.17 

15486.10 

8440.52 

7.850 

16275.3 

16275.3 

22918.7 

60 

2310.65 

16251.14 

7675.48 

7.850 

15010.8 

15010.8 

21262.7 

65 

2294.14 

17010.73 

6915.89 

7.850 

13712.4 

13712.4 

19532.7 

70 

2277.62 

17764.87 

6161.75 

7.850 

12380.5 

12380.5 

17729.5 

75 

2261.11 

18513.56 

5413.06 

7.850 

11015.4 

11015.4 

15853.7 

80 

2247.56 

19256.12 

4670.51 

7.850 

9618.9 

9618.9 

13908.0 

85 

2243.14 

19995.82 

3930.80 

7.850 

8186.2 

8186.2 

11886.7 

*70 

2248.46 

20735.90 

3190.72 

7.850 

6712.1 

6712.1 

97*3.1 

95 

2260.42 

21478.62 

2448.00 

7.850 

5192.4 

5192.4 

7591.2 

100 

2279.00 

22225.26 

1701.36 

7.850 

3624.2 

3624.2 

5306.9 

MAX 

2666.59 

22225.26 

17370.08 

7.850 

28072.6 

28072.6 

36323.3 

MIN 

2243.14 

6556.54 

1701.36 

7.850 

3624.2 

3624.2 

5306.9 

AVG 

2376.86 

14604.72 

9321.90 

7.850 

16892.4 

16892.4 

23147.1 

100* 

2276.10 

105 

1350.28 

22743.75 

1182.87 

2535.6 

2535.6 

3721.3 

110 

1052.03 

23136.73 

789.89 

1706.1 

1706.1 

2510.8 

115 

803.97 

23438.12 

488.51 

1063.3 

1363.3 

1569.1 

120 

584.86 

23666.32 

260.30 

571.3 

571.3 

845.6 

125 

375.92 

23821.76 

104.86 

232.1 

232.1 

344.5 

130 

157.53 

23908.10 

18.52 

3.3 

0.5 

0.7 

PER  CENT  PROPELLANT  (  0- 

100  PCT  WEB) 

BY  VOLUME 

»  65.49 

PER  CENT  PROPELLANT  1100- 

130  PCT  NEB) 

BY  VOLUME 

-  7.11 

TABLE  XII.  FINAL  PAGE  OF  OUTPUT,  TEST  CASE  1 
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THIOKOL-ELKTON 

_ ADVAiiCiD  MAIN  DESIGN  COMPUTER  PROGRAM . 

St  raight  Thru  Section  Broken  Back  Star  Variable  Temp  Option  Employed  to  Read 

Table  of  Const  Grain  Temp 

THEORETICAL  ANALYSIS  BY  F.E. MOORE 
PROGRAMMING  ANALYSIS  BY  D.H. FREDERICK 


InEE 

BURNING 

CHAMBER 

PROP. 

PROP. 

PROPELLANT 

INERTIA 

TENSOR 

rtURN 

SURFACE 

VOLUME 

VOLUME 

Cti 

I-XX 

I-YY 

I-ZZ 

PCT 

( SC— IN) 

(CU-IN) 

(CU-IN) 

(IN) 

- - - (  SQ- 

FT*LU— IN) 

0 

2634.05 

6556.54 

17370.08 

7.850 

28072.6 

28072.6 

36323.3 

5 

2666.59 

7430.92 

16495.70 

7.850 

27138.3 

27138.3 

35452.4 

10 

2620.76 

8310.01 

15616.61 

7.850 

26163.4 

26163.4 

34505.8 

15 

2530.06 

9159.67 

14766.95 

7.850 

25173.5 

25173.5 

33495.6 

20 

2445.54 

9978.20 

13948.42 

7.850 

24171.6 

24171.6 

32425.8 

25 

2A26..25 

10781*46 

13145.16 

7*850 

2 3145. Q _ 

23145* l 

.31289,4 

30 

2409.74 

11579.18 

12347.44 

7.850 

22085.2 

22085.2 

30080.0 

35 

2393.22 

12371.46 

11555.16 

7.050 

20991.5 

20991.5 

28796.8 

40 

2376.71 

13158.30 

10768.33 

7.850 

19863.8 

19863.8 

27439.3 

45 

2360. 2C 

13939.68 

9986.94 

7.850 

18701.9 

18701.9 

26007.2 

50 

2343.68 

14715.61 

9211.01 

7.850 

17505.7 

17505.7 

24500.3 

-55 _ 

_ 2327, 17 

.  15486.10 

..  944.0,52 

_7_.8S.Q_ 

16275.3  _ 

16275.3 

22918.7 

60 

2310.65 

16251.14 

7675.48 

7.850 

15010.8 

15010.8 

21262.7 

65 

2294.14 

17010.73 

6915.89 

7.850 

13712.4 

13712.4 

19532.7 

TO  , 

2277.62 

17764.87 

6161.75 

7.850 

12360.5 

12380.5 

17729,5 

75 

2261.11 

18613.56 

5413.06 

7.850 

11015.4 

11015.4 

15853.7 

80 

2247.56 

19256. 12 

4670.51 

7.850 

9618.9 

9618.9 

13908.0 

B5 _ 

2243.14 

19.995*82 

. 393Q.0Q. 

_ 7. 95.0.  _ 

8186.2 

0186,2 

_J1886.7 

9C 

2248.46 

20735.90 

3190.72 

7.850 

6712.1 

6712.1 

9783. 1 

95 

2260.42 

21478.62 

2448.00 

7.850 

5192.4 

5192.4 

7591.2 

100 

2279.00 

22225.26 

1701.36 

7.850 

3624.2 

3624.2 

5306.9 

MAX 

2666.59 

22225.26 

17370.08 

7.850 

28072.6 

28072.6 

36323.3 

HIM  .... 

_  2243. 14 

_ 6556  *54. 

..  J2.C1.35 

7,950 

3624,2 

3624.2 

5306.9 

AVG 

2378.86 

14604.72 

9321.90 

7.850 

16892.4 

16892.4 

23147.1 

100+ 

2276.10 

105 

1350.28 

22743.75 

1182.87 

2535.6 

2535.6 

3721.3 

110 

1052.03 

23136.73 

789.89 

1706.1 

1706.1 

2510.8 

115 

803.97 

23438.12 

488*51— 

_  1063, 3.  . 

1963,3 

1569.1 

120 

584.86 

23666.32 

260.30 

571.3 

571.3 

845.6 

125 

375.92 

23821.76 

104.86 

232.1 

232.1 

344.5 

130 

157.53 

23908.10 

18.52 

3.3 

0.5 

0.7 

PER  CENT  PROPELLANT  (  0-100  PCT  WEB)  BY  VOLUME  »  65.49 
PER  CENT  PROPELLANT  UOO- 1 30  PCT  WEB F BY  VOLUME  *  7.11 


TABLE  XIII.  FINAL  PAGE  OF  OUTPUT,  TEST  CASE  2 
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As  a  further  check  on  the  accuracy  of  the  program,  a  graphical  solution 
of  these  test  cases  (uniform  grain  temperature)  was  performed.  Figure  12  illustrates 
a  comparison  of  these  solutions. 

The  variable  temperature  test  case  used  the  temperature  profile  shown 
in  Figure  13.  This  motor  was  evaluated  by  both  graphical  and  computer  techniques 
and  the  results  are  also  shown  in  Figure  12.  Table  XIV  is  a  reproduction  of  the 
computer  output. 

The  variation  between  the  computer-derived  and  graphically-derived 
solutions  is  greater  for  the  variable  temperature  test  case  than  the  uniform  temperature 
test  case  as  expected.  The  principal  reason  for  this  is  that  any  small  errors 
occurring  in  the  graphical  analysis  of  the  variable  temperature  test  case  are 
cumulative,  while  those  occurring  in  the  constant  temperature  test  case  are  not 
cumulative. 

The  conclusion  of  maximum  importance  to  be  reached  by  comparison  of 
these  analyses  is  that  both  graphical  and  computer  solutions  follow  the  same  general 
trend  and  hence  tend  to  substantiate  one  another.  Further,  the  departure  of  the 
two  derived  results  is  relatively  small  and  well  within  the  errors  that  are  anticipated. 
The  effects  of  a  two-dimensional  variation  in  grain  temperature  may  now  be  analyzed 


by  the  Advanced  Grain  Design  Computer  Program. 


A4867 


FIGURE  12.  COMPARISON  OF  SOLUTIONS,  GRAIN  TEMPERATURE  TEST  CASES 
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THIOKOL-ELKTON 

_  _ _ -ADVANCED  DRAIN  .DESIGN  COMPUTER  PROGRAM . . . 

Straight  Thru  Section  Broken  Back  Star  Variable  Temp  Option  Employed  With 
Grain  Temp  Shown  in  Figure  13 
THEORETICAL  ANALYSIS  BY  F.E. MOORE 
PROGRAMMING  ANALYSIS  BY  D.H. FREDERICK 


WEB 

BURNING 

CHAMBER 

PRLP. 

PROP. 

PROPELLANT 

INERTIA 

TENSOR 

BURN 

SURFACE 

VOLUME 

VOLUME 

CG 

I-XX 

1-YY 

I-IX 

PCT 

(SO-IN) 

ICU-IN) 

(CU-1N) 

C  IN) 

- i  SO- 

FT*CU-IN) 

0 

2634.05 

6555.75 

17370.87 

7.850 

27254.5 

27254.5 

34686.2 

5 

2690.68 

7314.03 

16612.59 

7.850 

26418.0 

26418.0 

33878.5 

10 

2728.19 

0090.19 

15836.43 

7.850 

25529.7 

25529.7 

32987.5 

15 

2737.8 8 

8864. 8C 

15061.82 

7.850 

24603.9 

24603.9 

32019.9 

20 

2698.69 

9646.74 

14279.88 

7.850 

23625.4 

23625.4 

30955.3 

Z 5.  _ 

—  -2X11.24 

-3.0437.38 

134R9.24 

....  7.850  . 

22590.5  .... 

.  22590.  5- 

_  -25711.6 

30 

2794.84 

11239.82 

12686.80 

7.850 

21500.9 

21500.9 

28524.2 

35 

2856.82 

12084.50 

11842.12 

7.850 

20322.6 

20322.6 

27131.5 

40 

2904.70 

12976.34 

10950.28 

7.850 

19036.4 

19036.4 

25576.8 

45 

2945.73 

13914.81 

10011.82 

7.850 

17641.5 

17641.5 

23857.9 

50 

2995.13 

14895.57 

9031.05 

7.850 

16138.2 

16138.2 

21970.5 

__  -55  ... 

_ 1034.35  . 

-15930.31 

7996.25. 

. 7.85 Q 

14503.6 

14503.6 

19882.3 

60 

3078.53 

17011.65 

6914.97 

7.850 

12740.9 

12740.9 

17590.8 

65 

JG57.88 

18150.32 

5776.30 

7.850 

10833.7 

10833.7 

15075.8 

70 

2812.41 

19301.47 

4625.15 

7.850 

8832.5 

8832.5 

12387.1 

75 

2557.47 

20403.53 

3523.09 

7.850 

6822.3 

6822.3 

9624.3 

80 

2298.65 

21461.11 

2465.51 

7.85C 

4812.5 

4812.5 

6811.6 

_ as 

_ 1517.49 - 

22270.75 

1655.86 

7.85a 

....  3257.3 

325 7i  3 

_4625jl0 

90 

1151.91 

22844.63 

1082.00 

7.050 

2147.5 

2147.5 

3060.2 

95 

854.19 

23278.67 

647.95 

7.850 

1298.0 

1298.0 

1856.6 

IOC 

624.34 

23595.26 

331.36 

7.850 

669.9 

669.9 

961.7 

MAX 

3078.53 

23595.26 

17370.87 

7.850 

27254.5 

27254.5 

34686.2 

MIN  . 

624.34 

6555.75  . 

331.36 

7.850 

669.9 

.669*9. 

_  961.7 

AVG 

2461.20 

15250.84 

8675.77 

7.850 

14789.5 

14789.5 

19678.6 

100+ 

606.85 

105 

377.39 

23811.58 

115.04 

234.5 

234.5 

337.8 

11C 

103.60 

23916.80 

9.82 

20.1 

20.1 

29.1 

115 

n.nn 

23926.6? 

.  O.QC  . 

. .  Q.C- 

-O.Q  _ 

0.0 

120 

0.00 

23926.62 

0.00 

0.0 

0.0 

0.0 

125 

0.00 

23926.62 

0.00 

0.0 

0.0 

0.0 

130 

0.00 

23926.62 

0.00 

0.0 

0.0 

0.0 

PER  CENT  PROPELLANT  (  0-100  PCT  WEB)  BY  VOLUME  »  71.22 
PERCENT  PROPELLANT  (100-130  PCT  WEB)  BY  VOLUME  -  1.38 


TABLE  XIV.  COMPUTER  OUTPUT  FOR  VARIABLE  TEMPERATURE  TEST  CASE 
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III.  EFFECT  OF  STRAIN  ON  BURNING  RATE 

Results  presented  In  the  previous  quarterly  progress  report1  indicated  that  the  de¬ 
termination  of  the  effects  of  volumetric  changes  (Poisson's  ratio)  on  linear  and  mass 
burning  rate  provided  a  promising  approach  to  correlating  strain  effects  on  burning. 
During  the  previous  quarter,  it  was  shown  that  strain  per  se  does  not  directly  influence 
linear  burning  rate  but  may  affect  other  material  properties  that  can  be  directly  related 
to  burning  rate  changes.  The  change  in  volume  with  strain,  also  referred  to  as 
"Poisson's  ratio,"  is  this  property.  During  this  quarter,  several  propellants  were  in-, 
vestigated  by  determining  the  effect  of  strain  on  Poisson's  ratio  and  then  relating  that 
parameter  to  any  changes  in  burning  rate.  On  this  basis,  a  criterion  for  this  type  of 
phenomenon  may  be  established  based  on  change  in  density  with  strain,  whereby  one 
could  predict  the  effect  of  strain  on  burning  rate  by  an  examination  of  Poisson' s  ratio 
at  that  strain . 

A.  Poisson's  Ratio 

When  a  material  is  stretched,  cross  sectional  area  changes  as  well  as  length. 
Poisson's  ratio,  V  ,  is  the  parameter  relating  these  changes  in  dimension  and  is  defined 
here  as: 

•y  _  Percent  change  in  width  _  A  W/Wp  ^  ^ 

Percent  change  in  length  A  L/Lo 

It  can  be  shown  that  if  the  volume  of  a  material  remains  constant  when  subjected  to  very 
small  strains,  y  is  a  constant  and  equals  0.50.  Generally,  materials  increase  in  volume 
when  subjected  to  a  tensile  strain  resulting  in  a  V  less  than  0.50.  For  most  materials 
V  lies  between  0.2  and  0  5  and  approaches  0.50  for  rubbers  or  liquids. 
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When  the  material  is  subjected  to  a  strain  sufficiently  large  to  prohibit  the  use  of 


infinitesimal  strain  theory,  2/  is  no  longer  a  constant  but  a  function  of  the  strain: 


V  = 


1 


/  s1/2  1 

■  (i)  ]■ 


-1/2 


a  -  l 

where  2/  =  Poisson's  ratio 


(1.2) 


a  =  Principal  extension  ratio  (1  +  ei) 

e 1  =  Strain 

V  =  Volume 

V0  =  Original  volume 

If  the  material  is  incompressible,  i.e.,  V/V0  =  1,  equation  (1.2)  reduces  to: 


f  -1/2 ' 

I  1-  ( 1  +  e  i) 


(1.3) 


It  can  be  seen  from  equation  (1 . 2)  that  the  lower  the  value  of  7/ ,  the  greater  the 
volume  increase  and,  hence,  the  density  change  due  to  the  strain,  e  ^ .  Using  these 
three  equations,  three  propellants  have  been  investigated  to  date;  and  a  fourth  propellant 
will  be  studied  during  the  next  period . 

An  apparatus  for  the  determination  of  the  volumetric  changes  of  propellant  has  been 
developed  and  consists  of  a  point  light  source,  an  end  bonded  cylindrical  tensile  specimen2, 
a  10X  lens,  an  Instron  tensile  tester,  and  a  sheet  of  graph  paper.  This  is  shown 
schematically  in  Figure  14 .  A  specimen  of  known  diameter  is  placed  in  the  path  of  the 
point  light  source.  The  image  is  magnified  ten  times  by  means  of  the  lens  and  projected 
on  the  graph  paper.  The  highly  magnified  diameter  is  measured  and  the  volumetric 
changes  are  calculated  from  the  measured  lateral  and  longitudinal  strains. 
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FIGURE  14.  APPARATUS  FOR  MEASURING  CROSS  SECTIONAL 
AREA  OF  A  SPECIMEN  DURING  A  TENSILE  TEST 
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Three  propellants,  in  addition  to  the  TP-H-1050  already  studied,  have  been  selected 
for  these  tests:  TP-J-3000  (plastisol),  TP-G-3133  (polyurethane),  and  a  modified 
TP-H-1011(hydrocarbon).The  latter  was  formulated  to  give  a  JANAF  modulus  equal  to 
approximately  1000  psl,  so  as  to  obtain  a  more  rigid  material  for  comparative  purposes. 

Poisson's  ratio  values  as  a  function  of  strain  are  shown  in  Table  XV  and  Figure  15. 
From  the  figure,  it  can  be  seen  that  TP-H-1050  follows  very  closely  the  curve  for  an  in¬ 
compressible  material  obtained  using  equation  (1.3).  The  curves  for  the  TP-G-3133, 
TP-H-1011,  and  TP-J-3000  propellants,  respectively,  show  increasing  deviations  from 
incompressibility.  Poisson's  ratio  for  TP-J-3000  at  10  percent  strain  only  has  been 
obtained  to  date .  Additional  tests  will  be  made  during  the  final  quarter  so  that  this 
curve  can  be  more  fully  described.  The  data  on  TP-G-3133  and  the  modified  TP-H-1011 
represent  one  test  at  each  strain  level .  At  least  five  tests  at  each  strain  will  be  made 
on  each  propellant  before  the  final  curves  are  established. 

B.  Strand  Burning  Rates 

The  technique  employed  to  study  the  effects  of  strain  on  burning  rate  was 

3 

previously  reported.  Several  1/4  x  1/4  x  6  inch  strands  are  elongated  to  three 
different  strain  levels,  4,  8,  and  12  1/2  percent,  and  burned  at  nominal  pressures  of 
500,  900,  and  1200  psi.  The  strain  level  is  accomplished  by  using  a  constant  strain 
device  capable  of  independently  straining  two  strands .  The  specimen  ends  are  inserted 
through  Teflon  collars  and  bonded  so  as  to  be  similar  to  end-bonded  tensile  specimens . 
Once  the  desired  strain  level  is  achieved,  the  sample  is  coated  on  all  sides  with  Arm¬ 
strong  A-l  cement.  After  curing  the  coating  is  relatively  stiff  and  the  strain  level  in 
the  strand  is  maintained.  The  burning  test  is  accomplished  in  a  conventional  strand 


burner. 
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TABLE  XV 

POISSON'S  RATIO  VERSUS  STRAIN 
Test  Temperature:  77  *F 


Percent 


Strain 

TP-G-3133 

TP-H-1011 

TP-J-3000 

TP-H-1050 

2 

0.479 

0.478 

- 

0.496 

4 

.0.470 

0.445 

- 

0.485 

6 

0.465 

0.425 

- 

0.463 

8 

0.460 

0.416 

- 

0.459 

10 

0.4553 

0.409 

0.325 

0.455 

12 

0.451 

0.405 

- 

0.462 

14 

- 

- 

- 

0.464 

18 

_ 

_ 

0.456 
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FIGURE  15.  POISSON'S  RATIO  VERSUS  STRAIN 

Strain  Rate:  2 . 0  in. /min 
Test  Temp:  77* F 
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The  results  of  the  burning  rate  tests  are  given  in  Tables  XVI  and  XVII  and  illustrated 
in  Figure  16.  The  data  on  the  plastisol  propellant,  TP-J-3000,  were  obtained  from  a 
Thiokol-Redstone  Division  report^  and  represent  the  effect  of  approximately  10  percent 
strain  on  the  burning  rate.  The  latter  data  are  used  at  this  time  for  comparative 
purposes  only  since  it  is  planned  to  retest  this  material  at  this  facility  during  the  next 
quarter  in  accordance  with  the  above-mentioned  test  procedure.  The  effect  of  the  strain 
on  the  burning  rate  of  TP-J-3000  is  clearly  shown.  At  1000  psi  pressure  there  is  an 
increase  of  33  percent  or  about  3.3  percent  for  each  percent  strain. 

TP-G-3133,  whose  Poisson's  ratio  deviates  only  slightly  from  incompressible 
behavior,  also  exhibits  a  strain  sensitivity  on  its  burning  rate  but  to  a  much  lower 
degree  than  the  plastisol.  The  increase  in  burning  rate  for  this  formulation  at  a  pressure 
of  1000  psi  at  10  percent  strain  is  about  4.3  percent  or  0.43  percent  increase  for  each 
percent  strain. 

Burning  rate  tests  on  the  modified  TP-H-1011  have  not  yet  been  run.  However, 
since  its  Poisson's  ratio  at  10  percent  strain  falls  between  that  of  TP-G-3133  and  the 
TP-J-3000,  its  strain  sensitivity  on  the  burning  rate  should  fall  between  0.43  and  3.3 
percent.  These  results  will  be  obtained  during  the  forthcoming  quarter. 

C.  Effect  of  Strain  on  Mass  Burning  Rate 

The  effect  of  strain  on  mass  burning  rate  is  also  being  considered.  When  a 
strain  is  imposed  on  a  strand  it  will  produce  an  increase  in  the  time  necessary  to  com¬ 
pletely  burn  the  propellant  (i  e.,  a  10-percent  strain  would  cause  a  10-percent  longer 
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FIGURE  16.  BURNING  RATE  VERSUS  PRESSURE 
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burning  time  for  an  incompressible  material  exhibiting  no  linear  burning  rate  dependency 
on  strain) .  The  mass  burning  rate  is  thus  reduced  by  the  time  increase  required  to  con¬ 
sume  the  sample.  Therefore,  mass  burning  rate  is  a  direct  function  of  the  reciprocal 
of  the  strain  for  the  incompressible  case . 

The  linear  burning  rate  generated  from  TP-H-1050  has  indicated  that  this  linear 
rate  is  independent  of  strain  up  to  12  1/2  percent.  Therefore,  a  12  1/2-percent  strain 
on  a  strand  would  cause  a  12  1/2-percent  longer  burn  time. 

When  the  linear  burning  rate  changes  with  strain,  as  in  the  case  of  TP-G-3133,  a 
more  complex  relationship  exists.  Considering  the  Poisson's  ratio  values  of  the 
modified  TP-H-1011  and  TP-J-3000,  it  is  expected  that  a  similar  relationship  will  exist 
due  to  a  greater  sensitivty  of  the  burning  rate  to  strain.  These  propellants  will  be 
analyzed  as  soon  as  all  the  burning  rate  data  are  available . 
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IV.  PROPELLANT  SLUMP  ANALYSIS 

Principal  effort  in  this  area  during  the  past  quarter  has  been  directed  toward 
the  development  of  an  elastic  analysis  for  a  finite  length,  cylindrical-type  grain 
possessing  a  spherical  head  end  and  undergoing  an  acceleration  in  the  longitudinal 
direction  (see  Figure  17).  An  approximate  numerical  technique  originally  developed 

c  6  7 

by  Southwell  and  more  recently  extended  to  rocket  grains  by  Parr,  ’  Messner  and 
Shearly11  was  found  to  be  especially  applicable  to  the  solution  of  this  problem. 

The  pertinent  equations,  boundary  conditions,  and  finite  difference  approximations 
are  described  below. 

The  initial  step  in  any  stress-strain-displacement  analysis  of  a  solid  propellant 
rocket  grain  is  to  obtain  a  solution  of  the  elastic  field  equations  subject  to  given 
boundary  conditions  and  a  prescribed  external  loading.  In  particular,  if  the  grain  is 
symmetrical  with  respect  to  its  longitudinal  axis  and  is  under  a  constant  acceleration 
load  in  the  axial  direction,  the  governing  system  is  considerably  simplified. 

There  are  a  number  of  methods  that  may  be  used  to  solve  this  system®’  ®.  However, 
all  have  inherent  difficulties.  This  may  be  attributed  to  the  fact  that,  in  essence,  the 
mixed  boundary  conditions  typical  of  this  class  of  problems  must  be  satisfied. 

In  the  following  analysis,  it  will  be  assumed  that  the  plane  lamina  of  Figure  17  is 
an  axial  cross  section  of  a  cylindrical  grain  capped  by  a  spherical  head  end  and  having 
a  straight-through  port.  It  will  further  be  assumed  that  the  propellant  charge  is 
linear  elastic,  isotropic,  and  homogeneous,  and  is  bonded  to  a  rigid  case  on  its  outer 
periphery  (C')  and  head  end  (D').  The  stresses,  strains,  and  displacements  to  be 


Direction  of  Acceleration 


FIGURE  17.  CYLINDRICAL  GRAIN  WITH  SPHERICAL  HEAD  END 
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determined  are  those  induced  by  a  constant  acceleration  load  applied  in  the  longitudinal 
direction. 

A.  Stress  Functions 

It  has  been  shown5* 10  that  elastic  bodies  having  forms  that  are  solids 
of  revolution  and  loaded  so  that  axial  symmetry  is  maintained  in  the  resulting  deformation 
are  amenable  to  numerical  solution  using  stress  functions  and  relaxation  techniques. 

In  particular,  for  the  conditions  described  above,  if  two  functions,  $  and"^ ,  of 
the  spacial  coordinates  r  and  z  (but  not  If  )*  are  defined  in  such  a  way  that  they  satisfy 
the  differential  system 


az 


-  I  8  $  +  92  $  = 


8H 


9^£  -  - 
9  r^  t 


9  r 

aSL 

r  9  r 


9  zl 


(2.1) 


+  1 
9  z 


2tJ-  _  82 


9  zz 


then  the  equilibrium  and  compatibility  conditions  are  satisfied  if  the  stresses  are 
defined  by5 


°r  =  — 
r 


r  9 : 


"z  • 


L  *  ... 


r  9  r 


(2.2) 

(2.3) 

(2.4) 


and 


*  See  Nomenclature,  p.  65 
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£  =  1  8# 

rz  r  0z  (2.5) 

Iz  is  the  effective  body  force  due  to  acceleration  and  is  equal  to  the  load  in  "g"s 

multiplied  by  the  mass  density  of  the  propellant. 

The  associated  strains  and  displacements,  expressed  as  functions  of  £ 
and  their  partial  derivations,  may  be  obtained  from  the  axisymmetric  stress-strain 
and  strain-displacement  relations  of  elasticity.  For  reasons  of  brevity  these  are  not 
included  here. 

B.  Cylindrical  Grain  with  Spherical  Head  End 

Using  standard  notation6*7,  we  define  (see  Figure  17) 
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in  order  that  the  governing  equations  (2. 1),  stresses  (2.2)-(2.5),  strains  and 
displacements  may  be  expressed  in  a  more  compact  nondimensional  form.  These 
may  now  be  written  as 
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The  simultaneous  solution  of  (2.6)  subject  to  the  boundary  conditions 
described  above  and  illustrated  in  Figure  18  constitutes  a  description  of  the  stress- 
strain-displacement  state  of  the  propellant  charge.  The  boundary  conditions  are  delineated 


below. 


Side  A  (p  =  a  ] 


On  this  boundary  the  normal  radial  stress  and  the  shearing  stress  are 


identically  zero  and  hence 


ap  =  0 


(2.18) 


"Plj  =  0 


(2.19) 


Side  B  (»=  H) 


On  this  free  boundary  the  normal  axial  stress  and  the  shearing  stress  are 


identically  zero  and  hence 

=  o 

Tpt^  =  o 

Side  C  (p  =  1) 


(2.20) 

(2.21) 


On  this  fixed  boundary  the  radial  and  axial  displacements  are  identically 


zero  and  hence 


Up  =  0 


(2.22) 


0 


(2.23) 


Direction  of  Acceleration 


FIGURE  18.  DIME NSION LESS  COORDINATE  SYSTEM 
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SLde.Pfr  =.f  J 

On  this  curved  spherical  boundary  the  radial  and  axial  displacements  are 
again  zero  and  hence 

u p  =  0  (2.24) 

and 


u^  =  0  (2.25) 

Corner  Points 

At  these  points  the  boundary  conditions  for  each  pair  of  adjoining  boundaries 
must  be  satisfied. 

C.  Finite  Difference  Approximations 
1.  Governing  System 

Upon  considering  the  nodal  molecule  of  Figure  19,  it  is  readily  seen  that 
the  finite  difference  representation  of  (2 . 6)  in  relation  to  the  central  node  is 


h 

2"o 


(04  -  02) 


+  01  +  02  +  03  +  04  "  4  0o  " 


0 


W4-V2)  +  Vl  "Vi  +Y*  +  ^4  -4  t{>0-  01  -  03  +  2  0o=O  (3.1) 

The  question  arises  as  to  what  alterations  must  be  made  in  (3. 1)  when  one 
or  more  of  the  points  of  a  nodal  molecule  lies  outside  the  boundary  of  the  grain.  Such 
points,  aptly  termed  fictitious,  are  present  for  both  boundary  nodes  and  interior  nodes 
when  such  nodes  are  separated  from  adjacent  nodes  by  a  boundary. 

It  is  clear  that  if  the  curved  portion  of  the  boundary  lies  to  the  left  of  the 
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mid-section  as  it  does  in  Figure  18,  then  in  reference  to  nodes  centered  in  the  interior 
there  are  three  possibilities;  namely, 

(1)  Nodes  2  and  3  are  fictitious 

(2)  Node  2  is  fictitious 

(3)  Node  3  is  fictitious 

Consider  the  first  case,  which  is  illustrated  in  Figure  20a.  Without  loss  of 
generality  it  may  be  assumed  that  0  and  'f'  can  be  expanded  in  a  two-dimensional  Taylor 
series  of  the  form 


The  fictitious  nodes  2  and  3  must  be  replaced  by  the  boundary  nodes  A  and 
B.  If  it  is  assumed  that  the  origin  of  the  coordinate  system  has  for  the  moment  been 
translated  to  node  O,  then  the  points  1,  4,  A,  B  have  coordinates  (h,  0),  (0,  -h), 
(-6Ah,  0)  and  (0,  <5gh).  When  substituted  into  (3.2),  these  provide  four  linear 
equations  in  terms  of  the  first  and  second  derivatives  of  0  and  if/at  nodeO  as  follows: 
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XA  =  Xo-  <5Ah  (|2L\ 


+  llAhl  2 

2 


(3.2c) 


XB-Xo  +  pB>j(|f  )o  + 


(3. 2d) 


When  these  systems  are  solved  simultaneously  for  the  first  and  second  derivatives  of  0,  cr 
and  ^  at  node  O,  the  governing  equations  (3.1)  are  found  to  be  of  the  form 


Kx0!  +  K20B  +  K30A  +  «404  +  K50o  =  0  (3.  3) 

Ki  ^1  +  K2  f  B  +  K3  Ya  +  K4  ^4  "  K5  ~  "  *1*1  ‘  "  Klj  *c 

where  Kj  (i  =  1,5)  are  constants  to  be  computed  for  each  such  node.  They  are: 

ki  =  T^a 


Pq6  B<1+  6B) 


_  2p0  +  <5fih 


K A  =—~.SL 
4  P0(1+6  B) 

K5=  2 


1  ~  6 A  +  1  -  5B  +  _JL_ 
6A(1+6A)  6b(1  +i5B)  2Po 


Similar  expressions  for  case  2  (Figure  20b)  and  case  3  (Figure  20c)  may  be  derived  in  an 


analogous  manner. 
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2.  Stresses.  Strains  and  Displacements 

The  stresses,  strains,  and  radial  displacement  in  finite  difference  form 
(relative  to  an  arbitrary  node  O  were  found  to  be: 


(ffp)0  =2li0  [02  -04  +  V2 

'  Y4] 

(3.4) 

(  9^0  “  2h  P  0  t-  ^0  “  04  ]  + 

1 - 1 

O 

+ 

h-* 

1 

0 

(3.5) 

(CT  ^>0  =  ~2hP  o  ['t'  2- ] 

+  jj  z 
(° 

(3.6) 

^7'o-2hP0 

(3.7) 

(3.8) 

(6p)o=  jpo*^"  ^4)+  P0d-^)  (02-^4) 

(3.9) 

"  2h  [V0  +  (l- 

(ev)o  =  -  i^L[?2-Y4+*(02- 
(  2hPo 

h  )J  +  >j0z 

(3. 10) 

(UP  >0  =  ^  0  +  d  "  V  >0O]  -  2/P0  7oZ 

(3. 11) 

3.  Boundary  Conditions 

The  boundary  conditions  in  finite  difference  form  for  sides  A  and  B  are  the 
same  as  those  given  by  Parr”  for  the  free  boundaries  of  the  cylinders  he  studied.  They 


are: 
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Side  A  (  P  =  a  ) 

O 

II 

a-° 

(3.12) 

Va-  (‘Mi  -  JO-|.  [i+^j-]}  >o-0 

(3. 13) 

Side  B  (  =  H) 

yo=  (0 0  -  «2  )  HZ 

(3. 14) 

2  ^  3  +  2  C^4  "  ^2]  +  02  +  ^4  '  20O  "2  ^0  0 

(3.15) 

Side  C  (P  =  1) 

On  this  fixed  boundary  it  follows  immediately  using  (2.22)  and  (3. 11)  that 
one  condition  which  the  stress  functions  must  satisfy  is 

»o--rVy  &•-  Wf]  <31“» 

u  ^  is  also  identically  zero  on  this  boundary  and  since  is  invariant 
with  respect  to  7^  here,  this  implies  that  e  ^  is  also  zero.  Thus,  after  using  (3. 1)  to  eliminate 
the  fictitious  values  02,V  2  in  (3. 10),  the  second  condition  which  the  stress  functions  must 
satisfy  is  seen  to  be: 

y  +  (fj  +  2  -  (1  -J/)  0^  -  (1  -V)  (j)  3 

(h-3>  (3-17) 

+  2  V  04  +  2  (1  -  22J)  0O  -  h  ^0Z  =  ° 
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Side  D  [p  =  £(i^)] 

On  this  boundary  u  p  is  again  zero  and  the  first  condition  which  the  stress 
functions  must  satisfy  is  the  same  as  (3. 16).  Namely, 


0o 


_i_  [y0-^z_] 

\  -  V  1+2/ 


(3. 18) 


The  second  boundary  condition  (2 . 25)  when  applied  to  (2  ■  8)  does  not  present 
a  form  that  is  attractive  from  a  numerical  standpoint.  However,  since  uy  is  invariant 
with  respect  to  arc  length  on  this  boundary 

4^  =  °  (3.19) 


But 


8  u?7  =  8  u?  d  w 

8  s  (  9 ^  dt 

Hence,  using  (3.19)  and  (3.20) 


+  £u» 

9PC 


Ae_ 

ds 


(3.20) 


Now 


and 


8  u»  =  -  9  u 


X 


»„  =  .  =  -  Q-tafl  f  -UL  ♦  m  ±1}  .  »z 

fl  r>(  7  0  L  8  p  8  p  J  ( 


9  f 
9  u 


8  P 


7  -  fpy 


7 

9  up 


9 


f 


r.  ,  2  (luo  ±0L 

X  »  7 


8u 


7 


£_  =  (l±^)  r  _L^  +  (1-^)  ®j|-  2/p  Z 
p  L  8  >?  ^|J 


(3.21) 

(3.22) 

(3.23) 

(3.24) 

(3.25) 
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Noticing  that  on  this  boundary 

=  1  =JL 

d  ^  p  ( 


(3.26) 


and  employing  (3.22)  -  (3.26)  a  more  tractable  boundary  condition  is  obtained.  Namely, 

ML  +  y  i£  -  flUL  -  (1  _v)  M.J.  ilsJLL  -  -£_L_  T(i  -2/)  w+^1 

9p  dp  18^  ''  9^J  p  t  (l  +j/)  L  (  -I 

(3.27) 

Using  backward  differences,  it  is  found  that  in  finite  difference  form 
(4. 26)  may  be  expressed  as 


f0~V4  +  V  <0o-04>-  [W  *0~  *i}]‘  -1- 

Po 


p0hz 
(1  +*) 


0 


(3.28) 


Corner  Nodes: 

The  satisfaction  at  adjoining  boundaries  of  the  boundary  conditions 
defined  above  consists  mainly  of  algebraic  manipulations  of  said  relations.  For  purposes 
of  brevity  and  since  considerable  effort  is  yet  needed  in  defining  the  computer  logic  that 
must  be  developed  in  order  to  make  the  above  analysis  practicable,  these  conditions 
are  not  included  at  this  time. 


NOMENCLATURE 


r,  0,  z 

\ 

o,  e.  ^  - 

u 


Up 

W 

u! 

ar 


Cylindrical  coordinates 
Dimensionless  cylindrical  coordinates 
Radial  displacement  (in.) 

Radial  displacement  (dimensionless) 
Axial  displacement 
Axial  displacement  (dimensionless) 
Radial  stress  (psi) 

Radial  stress  (dimensionless) 

Axial  stress  (psi) 

Axial  stress  (dimensionless) 
Tangential  stress  (psi) 

Tangential  stress  (dimensionless) 
Shear  stress  (psi) 

Shear  stress  (dimensionless) 

Radial  strain  (in. /in.) 

Radial  strain  (dimensionless) 

Axial  strain  (in. /in.) 

Axial  strain  (dimensionless) 

Tangential  strain  (in .  /in . ) 

Tangential  strain  (dimensionless) 
Shear  strain  (in. /in.) 

Shear  strain  (dimensionless) 
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NOMENCLATURE  (Cont'di 

Body  force,  mass  density  x  acceleration  in  "g" 
Stress  functions 
Stress  functions 
Length  of  cylinder  (in. ) 

Case  radius  (in. ) 

Port  radius  (in. ) 

Elastic  modulus  (psi) 

Poisson's  ratio 

a 

b 

b  +  1 
b 

Arc  length 


Nodal  dimension 
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V.  FUTURE  WORK 

The  remaining  two-dimensional  heat  transfer  test  cases  will  be  completed  during 
the  final  quarterly  period.  This  will  effectively  complete  the  demonstration  of  this 
computer  program . 

The  capability  of  automatically  performing  a  complete  grain  design  analysis  with  an 
irregular  temperature  profile  will  be  installed  into  the  Advanced  Grain  Design  Computer 
Program  after  the  program  is  fully  checked  out  for  the  half  star  point  case. 

The  relationship  between  Poisson's  ratio  behavior  and  the  effect  of  strain  on  burning 
rate  will  be  defined  more  conclusively.  The  remaining  laboratory  testing  required  to 
obtain  a  significant  correlation  will  be  accomplished.  The  analysis  of  these  data  will 
consider  both  linear  and  mass  burning  rates. 

Immediate  effort  in  the  propellant  slump  area  will  be  directed  toward  the  definition, 
by  means  of  flow  charts,  etc. ,  of  the  computer  logic  necessary  to  obtain  numerical 
results  from  the  analysis  described  in  Section  IV  of  this  report. 

Complete  compilation  and  documentation  of  each  study  area  will  be  accomplished 
preparatory  to  the  program  summary  report  publication. 
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